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ABSTRACT 


Time~f>'(?'quc.ncy represent"' tion of siijnal and 

sysl-Gms i-s a poweilul tool in &lgnal processing specially 

when UiQ phenomenon unaor consideration is a nonstataonary 

one Such a reor’^sontation has many advent ages It simpli-* 

fie*- the analysis, helps in synthesis and design and gives 

a better insight about the signals, systems and their 

interaction Many such rei^re sen tat ions have been proposed 

Jn the p ist and have been found useful in certain aoplica- 

tions Recontly one such distiibution 'Wignsr Distribution’ 

(/©) has abtracteu the minds of scientists from various 

firlds like mathematics/ optics and signal processing mainly 

because it has many interesting properties This thesis tries 

to study the WD as a member of Cohen's class of time-frciiuency 

distributions. Mathematical properties of the WD have been 

listed and its applications to the fourier optics have been 

discus sod in brief A nev^ scheme of reception of pulsed 

chirp signal based on the VJP has been proposed and its per for- 

oi 

mane© i*i compared with/the correlator receiver. The now 
receiver performs better in detecting the presence of a weak 
return falling in the side-lobe of a strong return Ho;/over 
the cost of this improvoment is paid in teOTS of the deyrada-. 
tion of output SNR 
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CHAPTER 1 


INTRODUCTION 


in many cases of practical Importance mere time 
domain (sign il) or frequency domain (spectrum) representations 
are not adequate to visualize the phenomenon under invest! 
gation and scientists are interested in having a mixed tlmo- 
frcqucncy representation of signals and systems* Spectrograms 
used extensively in speech-analysis and the Ambiguity lunation 
used xn Radar are some such exaitples* Such a representation 
his many advantages It simplifies the analysis of the system# 
helps in synthesis and design and also gives a better insight 
about the signals, systems and their interaction. This is 
particularly so if the phenomenon under consideration is a 
nonstationary signal which means that at different instants of 
time the signal has different spectral oontent. 

Through the Fourier Transform, frequency is defined 
only over an infinite time interval and time is defined only 
over an infinite frequency band Thus it would appear, at 
least apoarantly /meaningless to ascribe the signal a combined 
behaviour. But this is not so Each 'line* of the frequency 
spectrum represents a continuous wave signal of that frequency 
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and phase and it 1 s to be sjqpected that at times i group of 
such continuous wave signols will interfere constructively 
and cause a bunching of signal energy within the spectral 
baiid of tlie group Thus initutively it does appear that 
energy muw>t have distribution in both time and frequency 
aom)Din( d 


Many suen time- frequency representations have been 
studied and Uued in the past and have boon found to be useful 
in certain ap£)lications All such representations however 
can be unified in a unique way Recently one such representa- 
tion the’ Wigner Distribution' (WD) has been tlie cynosure of 
scientists from different fields such as oiitics/ signal 
processing and mathematics mainly because it exihibits many 
interesting properties The concept of the V/D was introduced 
in 1932 by the famous physicst S Wigner for describing the 
probability densities in quantum mechanics fl] It was 
revived by Ville in 194-8 [2] and was used for the first time 
in signal processing but was again side tracked till it was 
found to be extremely useful in optics recently [3~7 3 

For dealing with stationary signals it is possible 
to estimate the energy content of the signal at different 
frequencies by taking its Fourier Transform and squaring its 
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modulus Value But this technique fails when the signal is 

nonstationary Tor such cl iss of signals what one needs is a 

representation which reveals how the energy is distributed in 
time- frequency combined domain An interesting example can 
be given which brings forth the limitations of the spectral 
density technique of analysis [8] Gon-^ider two signals 

S 2 ^(t) =Tl^(t-T/2) sin f 2jrt (f^-aP+at)] 

_ /.V sinTiBt . 

So (t) « QQg 2T[f t 

^ G 

where B » 2a T defines the bandwidth and the central 

c 

frequency, "• §) Indicates rectangular pulse with width T 
centred at t/ 2 rig 1,1 shows their spectral densities which 
are quite similar This is because the phase information has 
been lost in this representation. But the fact is that Sj^(t) 
is linearly frequency modulated signal and Sp(t) is not. This 
clearly shows the need for a tool which will give a better 
insight into the signal. 

The Time-rrequency Distribution (TrD) approach helps 
in solving such problems Many definitions of TFDs have been 
proposed with different kernel functions all of which have 



c 



c> 
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similar mathematical basis although the details differ This 
thesis tries to study a general class of Time-t requency 
Distributions vis a vis detailed study of tho WD The outline 
of the thesis is as follows 

Chapter 2 gives a general class of Time-rrequency 
Distributions with specific examples and their coirparisons 
Chapter 3 deals with mathematical properties of the WD 
Chapter 4 mentions some of the application areas of the WD in. 
optics Chapter 5 introduces a new scheme of detection of 
pulsed ohirp signals based on the WD and finally chapter 6 
includes suiwnary, conclusions and remarks 



CHAPTER 2 


GENERAL CLASS OP TIME PRCQUENCY DISTRIBUTIONS 

2 1 INTRODUCT ION 

As has been mention od in chapter a mixed time 
frequency representation or distribution has mmy advantages in 
unuerstandiny the behaviour o£ '’ignals and systems as well as 
in their analysis Many such distributions have been suggested 
in the past [9] These distributions generally have different 
properties \n efficient way to oonpare them systematically 
is to cozisider a general class of time- frequency distributions 
that include them all. Section 2 2 gives one such class 
called Cohen’s class of time*. frequency distributions Section 
2 3 mentions some members of this class and also oonpares 
them with each other 

2 2 CO HEN’S CL ASS OP TI^ -PREQUC NGY^ DlBTRIffl TlpNS 

A general class of time-frequency distribitions was 
introduced by Cohen [10], [11] a^^d was later widely rewritten 
by p Plandrin and B Esoudie [12] >[15] * They showed that all 
the possible definitians of the time-frequency distributions 
may be unified in a general definition 
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4-00 4 -^^ 

^ ^ 5 ^u) 

^OQ ^0® ^oe 

f (u-i^) f* dudTdS 

where f(fc) is the time signal# f (t) is its complex congugate 
and (5 #t) is a 'kernel function# representative of the 
partiaular distribution 

We wish to give a particxilar distribution an 
interpretation of energy density over time and frequency 
This imposes certain constraints on the kernel function 
®( S / t) Some other constraints are required to make the 
distribution meaningful in signal analysis For the puroose 
of systematic study and conp arisen of all the distributions 
belonging to Cohen's class a suitable set of properties was 
proposed by Classen and Mecklenbraukor [14] Those properties 
and the corresponding constraints on the kernel function are 
listed in Table 1 [15]* 

The first two properties are very useful because 
they enable us to give the distribution an interpretation of 
energy density over time and frequency. This is because 



Table 2.1 


Dxflerent properties Pj^and the coi re spending constraints 
on the kennels P(w) is the Fourier transform of the time 

signal f(t) 



Properties 


Constraint 

0 Kernel 

pi 

kn • 

= lf(t)p 

0) 

Kj ; 

1 fox 

• all 1 


4-00 

1 » 0 






P2 

J" C^(t, W )dt « 1 

^OO 

\rw\ 

0(0, X) 

= 1 

for 

all 

t 

P3 

If g(t) =s f(t ~ 

t.) 

00 . -c) 

does not depend on t 

p4 

then C (t, « ) = 

If g(t) ^ f(t)e ° 

0(^ /t:> 

does not deoend on m 


then C^(t,w)aC^ 

(t/ 0 ) - 0 )^) 






P5 

CfCt/ 

0 )) 

0(^ tX) 

sc ~ fc j 

T ) 





+ 00 


#'t)dS 


for| T 1<2 jt| 

P 6 

If f(t) =• 0 for 

{t|>T 

/e 

0(1 

« 0 


then C^(t, O') « 

0 for }t|>T 





P7 

If r(05) sj 0 for 

}o>l >Q 

/ e 

0 ( 

^ , i:)dT = 

= 0 for|J3-j<2|w| 


Then G^(t, w) » 

0 for| 0 )j >0 

mOO 






l-co 








/iC£{t/ w) dt 


0 ( 0 , C) 

« 1 

for 

all 


P 8 

_oo 

-Poo 

/ C^(t,0) ) dt 
^00 

=» t ^ 0 )) 
y 

— ^«((S 


^ =0 

S3 0 

1 

for all T 


•400 


0 (^/ 0 ) 






SUiC^itf w) dO) 


= 1 

for 

all 

£ 

P9 

*-00 

4.C0 

/'C.(t,w )dW 

= Q (t) 


x) 


cs 0 

for all ^ 1 


^00 



T =0 




plO C^(t, w )^0 for all t and is the ambiguity func 

^ tion of some functionh(t) 
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|f(t) I ^ 


is called energy density at time t since 


(a 

S |f(t) I ^ dt » Total Energy 

CO 


2 

Similarly |P(^>| is callv-d energy spectral density at 

u) 

frequency ^ since 

A 

\r po 

/ jF(<^')|^ dO) = Total Energy 

00 


Thus C^(t, w) admits of an interpretation as energy 
density at time t and frequency w if the integration of C£(t,w) 
over all frequencies at a fixed time t is the energy density 
at that time and the integration over all times at a fixed 
frequency w is the energy spectral density at that frequency, 
This in turn guarantees tnat 

,^ 0 O *^00 

1 

^ / / C^(t, (i)) - Total enexgy* 

^00 ^ CO 

as required. 

Properties P3 and P4 state that shifts in time and 
frequonoy give corresponding shifts in the distribution The 
next property P5^ which is very convenient from a practical 
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point of viBff, is that the distribution is real valued. The 
finite support properties P6 and P7 are important Th^ state 
that if a signal is bounded in time or frequency then the 
distribution will also be bounded in the same tune or frequency 
band 


The next two properties can be very useful for 

signal analysis Property P8 has the consequence that the 

^ . clUTiq ho'i'i 

centre ox gravity or average in the time 'direction at a fixed 

frequency of the distribution of the impulse response of a 

linear time-invariant system is equal to the group delay of 

the system at that frec[uenay The group delay is defined as 

usual and is equal to 

where ^2^(w) is the phase 

Property P9 is a converse statement. 

The last property PlO guarantees the positivity of 
the distribution for all times and frequencies For a meaning 
ful interpretation of o) ) as energy density at time t 

and frequency w this property is a must 

There has been a lot of discussion regarding the 
marginal ity requirement (properties Pl and P2) and the positivity 
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requirement (Property Plo) of the distributions Wigner has 
shown that for bilinear distributions it is not possible to 
satisfy both the requirements [16} This is due to the cons 
traint of Heisenberg's uncertainty principle Although 
Wigner was quite clear on the bilinearity condition, many 
subsequent authors have assumed that positive distributions 
satisfying marginality conditions do not e tist at alX, This 
has led to atbempts to give some interpretation to the negative 
values of the distributions and has introduced artificial 
and false concepts Nonbilinear distributions satisfying 
both the constraints do exist and have been put forth by many 
authors [X7J Recently Cohen and Posch have given an extremely 
useful and generalized procedure of generating non-bilinear 
distributions satisfying both the requirements [it] This 
in no-way contradicts the Heisenberg's uncertainty principle 
because the proof of the principle is based only on the margi 
nality property and on the fact that the distribution is 
bilinear [19] 

However till today mast of the distributions used 
commonly in signal processing are bilinear ones and hence can 
satisfy only one of these two requirements Spectrogram , 
defined in the next section and used widely in speech analy 
sis [20]^ is ah example which satis positivity condition 
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but not the marginal ity conditoLon In fact it has been shown 
Lhab spectrogram is the only example of Cohen's class with 
positive distribution [2l] Ambiguity function/ Jigner 
Distribution/ Levin Distribution^ all defined in section 2 3/ 
satisfy marginality condition but not the positivity condition 
The choico whether a distribution to be used should have the 
positivity property or the marginality property, depends upon 
the application in hand But generally the marginality 
requirement dominates in most of the distributions used in 
signal processing and hence negative values of the distribution 
must be accepted and interpreted. The most common interpretation 
of the negative values still assumes the distribution value as 
the energy Penalty at time t and frequency w in the sense that 



gives the energy of the signal in the time interval t - 
to t 4 and in tho frequency band <0 ~ ^ to o) I- 

But time interval At and the frequency band Ay cannot both be 
made arbitrarily small simultanoously. Uncertainty principle 
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restricts the product /Su to a lower bound and it can 

bo shown that when At Aa> exceed this bound, the above 
double integration always yields a positive value [ 22] 

Hence C^(t, o ) rnay have negative values at certain points in 
the t~0) plance but the total energy over any small area of the 
t -Cl) plane is always positive, as long as the aroa used for 
corrpubation is greater than some minimum value as given by the 
undertainty principle This explaination is intitutively 
satisfying and removes the erribarassment oaused by the negabive 
values, 

2 • 3 SOME K NOWN TIMC -FHEp UENC Y PIS TRlBUriONS 

Some of the Joiown time-frequenay distributions with 
their kernel functions have been listed in Table 2 fl5] 

This table also includes the properties satisfied by a parti- 
cular distribution It can be seen from the table that many 
of these distributions satiofy all the properties except the 
positivity property Hence bo choose between them we ne'^d an 
additional criterion. An important criterion is the spread 
of tho square of the magnitude of ) [<25] As an 

example let us compare two distributions, namely Wigner 
Distribution and rc al part of the Rihaczek distribution, for 
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a chirp signal f(t) given by 


X5 


f(t) = 

Tne insl-antaneous frequency of this signal is given by 2at, 
The Wigner Distribution of this signal is equal to L15] 

W£(t, w) = 27T5(w-2at) 

This result is extremely satisfying becau^je it 
gives nonzero distribution value at time t only for frequency 
aa2c.t which is the instantaneous frequency of f(t) at that 
time 

The Rihaezek distribution can be shown to be equal to 

C£(t, w) - - j |) 

(the real part of the Rihaezek distribution in which we are 
interested is equal to 

CfCWw ) cos J) 
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This shows a large spread around the lino w = A better 

comparison ■’would be to comoare the spread S(t|^/ of the 

square magnitude of ) at a point (t^z of the t- w 

plane defined as 

y 20 4-00 

S^(t^z = / /((t-t^)^h(0>-W^)^) jC£(t,6)) . 

„ OO „ 00 

It Can be shown that the spread of the square magnitude of the 
Wigner Distribution is minimal as conpared to the spread of the 
Rihaczek Distribution! 2^] Hence Wigner Distribution is 
better of the two in analysis dealing with tho cnirp signals 



CHAPTtlR 3 


WIGNCR DISTRIBUIION INPRODUCTIOa AND PROPLRTIDS 

3*1 DEriNITION OF JHE WIGNER DI&IKIBU TION 

In, chapter 2 the defxriition o£ the Wigner Distri- 
bution (^D) as a member of Cohen's class of distributions/ 
has been gxven More generally the cross ^\ID of two time 
signals f(t) and g{t) is defined as 

0Q 

(t/ti)) « ; f(tr7?/2)g’''(t-r)/2)e"^‘^^ dH (3 1) 

# y 

--00 

where the ’’ ^ denotes conplex conjugation 
If f{t) = g{t) we get auto-WD of the signal f(t) 

bet P(w) and G(W) be the Fourier Transforms of f(t) and g(t) 
respectively A similar expression for the cross v® of the 
two spectra P(^) and G(W) is given by 

4 * CO 

Wj,^Q(W/t) F(4) +V2) G»^(W-V2) d^ (3 2) 

M 00 

It can be shown that 




(3 3) 
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This illustrates the symmetry between the time and frequency 
domain definitions of tho 

The term •Distribution’ in the ncime vVigener Distri- 
bution suggests the WD to be similar to a joint probability 
density for the ‘random variables’ time and frequency The 
reason for using probabilistic language (e g distribution, 
expected value, marginal ity etc ) in discuising the WD or any 
other time- frequency distribution for that matter, is two - 
fold There is a historical reason that these distributions 
were first used in quantum mechanics which is inherently a 
probabilistic theory Hence the original papers by Wigner, 
Ville or Cohen contain methods and language of quantum 
mechanics The theory of the time-frequency distributions as 
established by subsequent workers has the same mathematical 
structure as joint distributions in probability theory In 
particular, the loquiremenbs one imposes on the distribution 
function (marginality, positivity) and the methods one uses to 
calculate physical quantities, that is, integration over the 
time- frequency plane with a density function is identical to 
the method used in joint probability theory Hence one may 
think the time-frequency distribution as a probability distri- 
bution and consider time and frequency to be random variables. 
Obviously, this is not essential always but is suggestive of 
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the concept involved and hence very useful 

In chapter 2 some of the important prop erties of 
the WD have been listed in Table 2 Section 3 2 of this 
chapter gives an exhaustive list of the properties of the WD 
Sections 33, 34 and 3 5 deal with the transformations of the 
t-Q plane which are invariant to the WD Sections 3 3 and 
3 i describe such linear transformations whereas section 3 5 
describes nonlinear transformations Section 3 6 presents 
some illustr itions of the WD, 


3 2 PROPERTIES OP THE WIGNER DISTRIBUTION 

Before lifting the properties of the v7D the 
notation used in the discussion is introduced here 


The Fourier Transform of a time signal f(t) is 

denoted as 


h*> 

r(o>) - P T (f(t)) = / f(t) dt (3 4) 

The inverse Fourier Transform of the spectrum 
r(w) is denoted as 

H CO 

f(t) « F T “^(F(0)) - S F(W) (3 5) 

*-^00 
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The -time signal f(t) and its as defined In 
(3 1) with g(t) » f(t), is denoted as 


f(t) »<• — ^ ^(t/Oi) 


(3 6) 


The inner product of two time signals C(t) and 
g(t) is defined as 




<^,g^ ~ ^ f(t)g’<(t) dt 


(3 7) 


Similarly the inner product of two spectra P(w) and G((i)) ii 
defined as 


<F,G> 




/ r(o5) G*(0)) dw 


27t 


-»*00 


The norm of a time signal f(t) is defined as 


1 
N 2 


(|f|| = <f'^> 


(3 8) 


(3 9 ) 


and similarly for the spectrum F(^) 


Ill'll .<P,F>Z 


I 


(3.10) 
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Several other operations on signals will be 
introduced below 

the shift operation 

S^,f « £(t-x) (3 11) 

the (complex) modulation. 

f(t) (3 12) 

dif ferenti ition 

Df = j 

multiplication by <-he running variable 

]jf=*tf(t) (314) 

and reversal of the running variable 

Rf « f(-t) (3 15) 

Most of the properties given belov/ are extenbions of 
similar properties of the Ambiguity Function [50] These can be 
proved easily from the definition of the WD and hence the. 
proofs have been omitted Some of the proofs and related 
discussion however can be found in [51] /f52i[53]and [14] • 

1, For any two time signals f(t) and g(t), 

_(t/W) =5 w* -;(t#(»>) (3 16) 

t/g g/x 

Hence the auto /TO of any (real or complex) signal is always 
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real Moreover the WD of a real signal is an even fmetion 
of the frequency/ i e 

(3 17) 


W. 


(t,w) *= (t/ -W) 


2 Time Shift A time shift in both the signals f(t) and g(t) 
corresponds to the same time shift in the WD/ i e. 




(3 18) 


3 Frequency Shift Modulating both the signals f(t) and g(t) 
by results in the same frequency shift of the WD/ i e 


W. 


M 


■F M ^(t,w-.a) 

Mjjg' f/g 


(3 19) 


4 Combining property 2 and property 3 we get 

This result can be used to express the WD as 
W. ^(t,o> “ W M f s 


(3,20) 


, ( 3 , 21 ) 


Using the definition of the inner product and the operator 
R given in (37) and (3 15) we get 


Wf^g(0,0) « 2 <f, 


(3 22) 
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From the above two equations it then follows that 

= 2 ( 3 - 23 ) 

Thus the value of the WD at a certain time t and frequency 
w can be determined by the inner product of the shifted 
and modulated signal# the second of which has undergone ^ 
time reversal 


5 Addition of Two Signals The ® of two signals f(t) and 
g(t) is a bilinear functional of f and g which means that the 
WD of the sum of the two sianals is not simply the sum of 
the WD‘s of the individual signals It can be shown that 


W 


^1 ’■^2 * 


91^92 


(t#o)) « vJ, 


^1 * ‘•^1 


(t# « )tW, 


1^^2 


(t#W)+w, 


2'^1 


(t,W) 


^ (t#W) (3 24) 

^2 '^2 

As a particular case 

(W. „(t#0))) *(3 25) 

y 


6 The product of the VJD of two signals f(t) and g(t) with 
running variable t can be e3q)ressed as a sum of two WD*s 


(3 26) 
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where the operator L has been defined in (3 14) Similarly 
the multiplication of the WD by w can be ejqpressed as 


where the operator D is defined in (3 13) 


(3 27) 


7 According to the definition (3 l),thc WD can be treated as 
the spectrim (or the Fourier Transform) of the signal 

T TT 

f(t-^)g*(t^^) considered as a function of x with t as a 
fixed parameter Therefore the inverse Fourier Transform 
yielids 

•\o« 


which Can be written as 

^ Jy-^,0)) dW = f(t^)gMt2) (3,28) 


If tj^«t 2 »t we get 


1 

2n 


+09 

; 


W£^g(t, ijs) dw * f(t)g*(t) 


"'00 

and in particular, 

^ 60 

/ W£^£(t, to) dw =* 



(3 29) 


JL 

2n 


•^00 


4 


(3 30) 
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This proves the mcirginality- property of the WD rnentioned 
in the Table 2 of ohaptex 2 

If t]^ « t and t^ » 0 we get 

W£^^(t/2, w ) d« = f(t)g (0) (3 31) 

•-* CO 

This IS an interesting result which shows that f (t) can be 
completely recovered from its cross WD at time t/2, upto 
the constant factor (0) Same is true for g(t) A simi- 
lar discussion applies to the definition of the WD given in 
(3 2), Here the WD can be considered as the Inverse Fourier 
Transform of the spectrum P(D+-)G’' considered as a 

function of ^ withti) as fixed parameter Therefore the 
Fourier Transform yields 

-f. 00 

Wj,^Q(W,t) dt - F((i)'- |)G’^{oJ -V2) 

00 

using (3 3) ahd changing the varoblias we get 

-t ft \ / \ N 4- 0)^ HO) o 

^ „ r(o^)GA((0^) (3 32) 

00 

If =ii >2 we get 
+ 00 

^ w. it, W ) dt = F(a)) G*(0)) 

g 


-* 00 


(3,33) 
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and in particular 

S dt = jr(«)|^ (3 34) 

— m 

This is another marginal ity property of the WD 
If 0)^ » w and (*>2 *= 0 we get 

J W. ^(t/ 0.'/2) dt - P(0))G’ (0) (3 35) 

y 

•-CO 

This relation shows that the spectrum of f(t) can be coirple- 
tely recovered upto the constant G* (0) from its cross WD at 
frequency w/2* 


8 Prom equation (3*29) it follows that 


1 


2n 


*Hn» 

/ 

-•00 


} 4 M 

^ „(t,0)) dO)dt 

r^g 

-•00 


/ f (t)g* (t>dt = 

(3 36) 


Similarly from equation (3 30) we get 
*b -to *b 

2 n ^ dt (3 37) 


Thus the energy of the signal between the time intervals 
t^ to tj^ oan be found from its WD similar results exist 
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for the spectrums F(w) and G(w) 

9 Moyal’s Formula Moyal's formula for the integration of 
the product of two WD‘s is as follows 


1 

2jt 


^ CO 4. 00 

^ „ (t#o>) Vfp (t/w) dtdw 



(3 38) 


lO, Tune Limited Signals If f(t) and g(t) are restricted 
to finite time interval only and 


f(t) M g(t) = 0 for ^ • (3 3^) 

then the WD is restricted to the same interval/ i e, 

W. (t/ w) « 0 for t < tcft / ^or all w ^(3 40) 

11 Band Limited Signals If f (t) and g(t) are both bands- 
limited/ then the WD is limited to the same band, i e if 

p(W) « G(o)) =5 0 forWs^W<0) C3 41) 

then 

W. (t/«) = 0 for “ (3 42) 

r/ g a o 

12 The WD and the Ambiguity Function As has been mentioned 
in chapter 2, the Ambiguity Function (AF) of two signals 
f(t) and g(t) is defined as 



2B 


o) = / &! (3 43) 

<>00 

or equivalently 

.|eo 

Afg(T, 0) / r(oi^)G*(^p.>|*)e-^’^^d« (3 44) 

mJM 

It can be shown that the WD and the AP of two signals are 
related by a two-dimonsional Pourier Transform/ i o 

hco 

Af^g(T./ 0) “ ^ ; e-J ^ 


and 




Wf „(t/ 0>) 

y 


1 

2rr 


-JeO 
/ 

>• 00 xoo 


/ e"*-^ ^ ^ ^ ^ dtdO 

y 


(3 46) 


13 If f(t) and g(t) aie both even or bot7:i odd functions of 
time/ ±*e 


f.<t) » + f(~t) (3 47) 

and g(t) » 4^ q(-'t)/ then the AP and •che WD are the same up to 
scale factors/ i o* 


W£^g(t/0)) » + 2 A£^g(2t/ 2W) 


(3 48) 


14 




All the T±me~Precxuenoy Distributions of the Cohen’s class 
Can be obtained by a linear transformation of the WD f 
characteriaed by the kernel <P(t, which is related to 
by a two-dimensional Fourier Ira ns form/ i e. 


J^CO 

C^Ct/ (*>) - 2 k ^ ^ 


g(T/1? ) dTd'P 

(3 49) 


l-CO *]00 

wh^ro cp(t/W) == 2? ^ ai(y ,T)d'I'dT (3 50) 

^co 

15 The auto WDS of two signals f(t) and g(t) are identical iff 

f(t) = g(t) (3 51) 

where X is a real constant 

16 If w) 4-> f(t) 

then W£^ ^(at/ w /a)-<^-)^^|a I f(at) (3 52) 

where a is a constant This means that the time stretching 
of the signal f(t) affects both the axis of the t-w plane 

17 If ^(t/to )4-^£(t) 

then ^(t/ 0) - 2bt)<-^f(t) = f(t)e^ (3 53) 

where b is a real constant 

This means that if the WD has most of its energy along the 
t-axis then by chirping the energy concentration is moved 
to the line o) » 2b t. 
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18 If Wg^(t^w) ^ ^ f.(t) 

then W££(t+ 2 dW, 0))^->'^(t) = F T ''^(F(co )) 

(3 54) 

where d is a real constant 

This moans if the has most of its energy along the axis 
then hy frequency chirping the signal, the energy concentra- 
tion is moved to the line t = - 2dW . 

19 hot ^Ct, w ) ^ ^ £(p) 

Wg^g{t,a) )^> g(t) 

and S(t, w )- w. .(t, W ) + ^(t, 0) ) (3 55) 

E/X g,g 

then S(t^ w ) IS an auto WD iff 

f(t) = k^g(t) ^ (3 56) 

and hence iff 

<3 b7) 

i ? 

where is any constant and k 2 - jk^[ 

2L Invariance \ander Rotation of t- Plane The rotated WD 
i^ still a VO which means that if 
W£^£(t, w) .^f<t) 

then 

Wf w cos O H t sin 0# Woos 0 - t sin 0)'^->-g(t) 


(3 58) 
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where g(t) ^ e 


. tanO 
2 


Tlj2 


/ p( 0) ) e 


. (0 tanO 
j -•-- 2 *’ 


e' 


jwtsecO 


dw 


21 Let W£^£(t, )<”> f-(t) 

„(t, w ) g(t) 

y / y 

)<~^h(t) 

a) If h(t) ” f(t),g(t) 


(i) 


then w^j,(t, CO ) = ) * Wg^g(t, O) 


0 ) 

where indicate convolution w r t w 


b) If h(t)=u(t) * v(t) 

then w) = W^^^Ct, O) Wg^g(t, w) 

where J indicate convolution w r t t 

22 Let u(t)/ v(t)/ g(t), h(t) be time signals 

w 

1 f t. W ) Vf vj f tiO)) 

““^2? V,ug ' 


“ 2n '^uh/Uh 


w 

(t, (»)) * W {t, 0)) 

^ * ' vg,vg^ ' ' 


(3 59 ) 

(3 60 ) 

(3 61 ) 

( 3 , 62 ) 

(3 63 ) 
(3 64 ) 


52 


x: 

and a> ) « W (t, W) w , r.(t, w) 

u.fv,g»h' * •ug^'ug' * ' vh,vh' ‘ 


or =: w 


'uh.uh<t' ■* %,V9''" 


(3 65) 
(3 66) 


7 3 Lot ^ (t, 0))^>f(t) 

Wg^g(t, w)<->g(t) 


0) 


Ihen 0.) . ^ o) * (-2) e 


- j 2 w t 


* 


(F T (§)) (-2W ) 


(3 67) 


where (FT (2) v) (_2(o ) indicate Fourier Transform of the 
time function g”(t)/f^(t) at frequency -20) 

24 Realiaability Condition of the WD A distribution C(t/C0 ) 
is an auto WD of some time signal iff 

^ 0 )) dW * f(x)f^(y) (3 68) 


where f(t) is any time signal and in that case ) 

C(t, 0) ) = W£^£(t, w) , (3 69} 

Similarly C(t^w) is an auto WD of some spectrum ifff 

-1-co 

/ C(t, dt = F(x) (y) 


— (3 70) 
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for some spectrum r( w ) and in that case, 

C(t, w) _ W££(t, u) (3 71) 

3 3 LINE AR TRAlBEORHATI QN OF THE t-.^ PLANE 

The WD of a given time signal may have high density 
areas in certain parts of the t-w plane If for some reason 
we wish to shift this area of high energy concentration to 
some other area, transformation of the t-w plane is called for 
But we are interested in only those transrormations which are 
invariant to the WD In other words the transformation snould 
be such that the distribution with transformed coordinates is 
still a V© corresponding to some time signal Thus if the 
transformation is invariant to the we can obtain a new time 
signal with desired distribution of energy in t-w plane from the 
given signal and its distribution of energy In this and next 
section we will deal only witn linear transformations of the 
t-w plane In this section we will assume linear transformation 
of o> axis only Next section deals with general linear trans~ 
formation of the t-oi plane. In section 3,5 we will discuss 
nonlinear transformations All this analysis is an extension of 
a similar analysis of the Ambiguity Function [50], 

Let the transformation of the t-m plane he as follows 
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» t (3 72) 

W' = ko) (3 73) 

where t‘ and w* denote the axes of the transformed 
plane and k is a real constant 

Let us see under what conditions this transformation leads to 
the invariance of the aTO In other words we wish to find the 
condition on the value of k such that after the transformation 
of the t-O) plane the modified time-frequency distribution remains 
a WD coriesponding to some time signal 

w ) > f(t) 

and C(t, w) - W) =» kw), (3 74 ) 

Prom equation (3 70),G(t/W) is also an auto WD iff 

/C(t, ^ )e'-Jt{x-y) ^ 


for some spectrum G(w ) 

i e iff / W£ j_ G(x)G*(y) 


S f (tH'j)£«(t.S)e 


_1 

~j 2 


-jt(x-y) 


dridt -G(x)G"(y) 


.J90 
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Let t+^ = p, t^^^ssq 

t ^ T) « p~q 

C(t, w ) is an auto WD iff 

Vfo ■\v> -j H-C-SiXl (p-,q) -j (x-y) (p lq)/2 

/ /f(p) f^' (q) e e dpdq 

00 — 00 

= G{x)G»(y) 

1 e Iff s{ ^ iSjSi] r _ i^-]= G(x)0' (y) 


i e iff F [(k'*l)x h {k-.l)y ] P’*f[(}c~l)x + (kl-l)y] = G(2x)G (2y) 

for dll values of x and y 
Putting X - 0 we get 

^L(k-.l)y] F*[{lc+l)y] « G(0)Gf(2y) 


Putting =3 0 we get 


J'r^c+Dx] P f(k-l) c] = G(?x)G*(0) 


Hence we get 


G(2x) G*(2y) = F[(k-l)y] F*{(k+l)y] 

^ LQ? hlA.^ 1 ^ 2 

|G<o)l^ 


C3 
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But |g(o)| ^ == I r(o)| ^ 

C(t/ <>') is an auto WD iff 

P [ (k;+l)x (lc~l)y] 1’ ( (]c-l)x + (l<-'3)y ] 


I P(o) I ^ 


LUs=l^J 

(3 76) 


This shov/s that only tnose values of k which satisfy equation 
(3 76) will give rise to an invariant linear transformation of 
t-w plane (assuming t* = t) It can be seen from equation 
(3,76) that k » +1 always satisfy this condition/ a result as 
eixpected. 

Special Case 

Let r(w ) be a rational function which is the most important 
case in practice/ i e 

7t((0 - Wj^) 

F(w) 

j ^ 

Let us find conditions under which G( w ) is also a rational 
function since this also is the most commonly required condition 
Now from ec{uation (3 75) with w ax and \=y we get 

G(2u) G*<2\) a P[(ki-1)« +(k-l)?v.] F' [(k~l) «+(krl)X ] 


(3 78) 
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n ( 0) - 0^ ) 

LSt G ( ^ ) = " » — -»'•> - » '.r*' . 

7T ( (j) “ 0) T / 

1 

Thus for all w=s(i)^y2 ^ ^ ® equation (3 78) above is 

zero for all values of except at \ = ^\/? 

F [O^+l) 1- (k-.l)\]= O (3 79) 

or (kU)X]= 0 (3 80) 

for all \*s except X m 

If F [(K.+1) ^ny2 ] = 0 for all values of \ except at 

X- ^^2,/2 must have k=sl and hence 

If ?>» [ (k~l) {khl)\ ]« 0 for all values of X except at 

N = then we must have k =» -1 and henco 

Thus for G(0) ) to be rational/ only tv/o values of k, namely 
k=jrl/ are allowed whersos if no condition is a.mpo'='ed on 
kae +1 are t-wo of the possible allowed valuer of k 


3 4 GENERAL LINEAR TR l^ SgORMATIOH OF gHB IlIi kM 

In this section we consider a general linear trans 
formation of the plane» 
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Let ^(-t, W ) (t) 

W. ^ (t, 0) )<->•£ (t) 

wheire f(t) a.s e time signal and ^ ^ signal such that 

W- ^ (t, w ) = _(t, ko) ) . (3 81) 


where k is a real constant 

Let g(t) = f (at) 

where a and b are real constehts 

Then W„ ^(t/ w ) ~ ^ (at, ~(tA) ~2bt) ) 

g.g ® 


— (3 82) 


, (3 83) 


by using equations (3 52) and (3 53) 

2 

Let h(t) = F T [G(W)e-^^“ ] 

where d is a real constant 

^ [ a(t r2dO ) , ~ (co -2b(t f-2d W ) ) ] ♦ (3 85) 

by using equation (3 54) 


This means that if we have a linear transformation of the t-w 
plane of the type 
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t' 


a 

2 ad 


~ t * 


i 

.^2bk 

a 

(l~4db)k 

■» *' fill* 

a 


0> 


then W£^£(t',w‘) « ^ 


(3 86) 


But for the existence of ^ (t^ O) ) and hence in turn for the 
existence of w ) , we need ]< such that 

P[ (k+l)w+ (k-l)\] F'[ (k-1) w +<ki-l)\ ] 

2 


P(0) 


(3 87) 


as given hy equation (3 76) with x andjr 


Prom the above discussion it follows that if we consider a 
general linear transformation of the t-w plane 


t‘ 


I 

O' 



^22 




^2 ^ 


t 

i 

^22 

i 

0) 





(3 80) 


then by aonp arisen we get 
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a 




b 


2(det«:>) 


k =* (dat, D) 


(3 89) 


where (det, D) denotes the determinant of the transformation 
matrix D 


The transformation will be invariant to the WD if (det D)-<'k 
satisfies the condition mentioned in equation (3 87) However 
since k = +1 are valid values of k we can claim that any matrix 
of transformation D is valid as long as 

det D - + 1 (3 90) 


3 5 NONL IHaA R TRAI^FO’^M ATION OP THE t>^t^ PLy\NE 

In this section we v/ill consider any general trans- 
formation of the t-w plane For the sake o± simplicity of 
analysis we will assume that only w -axis has been tranformad* 
The result for t-axis transformation follows exactly the same 
pattern and hence will not be discussed here 

Let o)'=iu(0i>) be a general transformation of the O) -axis 

Let ^ ^ 

C(t/6> ) = u(w)) 


and 


(3 91) 
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liet us find the condition under which the transformation u(w) 
is invariant to the WDj i e we wish to find condition under 
which C(t, w ) is also an auto WD of stme signal 

4 eo 

Let T^(x,y) » I c(t, (3 92) 

OO 

= / f(ti"^) f* (t-.~) e ^-^^(x-y) 

•^00 00 


Let tilV2 » p t-T} /2 = q 

4«. ho (p-q) 

T^Cx^y) = ; / f(p)f*(q) e 


g-j (x-y)(E^) 
dpdq 


-jeo >0 ** JP r 

I ! f(p) e 


M. 00 ^00 


jqf- 2C^ + u(^^)] 
f' (q) e 

dpdq 


= F [u + 2^] I* [ vuc^fC) _ ^ ] 


By equation (3 70) C(t/ 0)) will be an auto WD and hence the 
transformation u(o ) will be invariant to the WD iff 
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r . 2^] P [ (3 93) 

« G(x> G*(y) 

for all X and y and for some spectrum G(w) 

Special Case Let F({*>) be rational as given by equation 

(3 7/) which is the mo.-t commonly found case in practice 
rrom the practical point of view we will be interested in 
having G(W) also a rational spectrum Let us find conditions 
under which G(w ) remains rational 

Let us for a moment assume that C(w ) is also rational 

n w - OJjjj 

G(w) » (3 94) 

1 ^ 

* <^) " 0 for 'VAexcept K « 

hX •-'X 

F [u.(-?’'2-) ] ■> 0 . (3 95) 

w H X , 

or — ) — 3“ 0 .(3 96) 

V* X except X ~ ^ 
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Equation (3 95) irnplies 


^ ) + 


= w 


w„ + X 

Let -“-2“^^ » ^ 


uW =o>^. .‘4 


= ( -<^) + ^ 


HQ. 


(3 97) 


whore is a constant given by 


K 


dm 


0 ) 


0) 


n 


Thus for all possible values of we get a transformation 

which transforms the w -axis as given by equation (3 97) 

Similarly equation (3 96) implies 

f('3? ) == K^j - W (3 98) 

and we get another set of valid transformations 

Hence if we wish to have rational G(<^) a large number 
of allowed transformations exist given by equations (3*97) and 
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(3 98) although tho nature of these transformations is the 
same for all It is important to note that all these trans- 
formations arc linear implying that no nonlinear transfor- 
mation of the w-axis exist if we wish that a rational spectrum 
F(w) should give rise to another rational spectrum G(w) after 
the transformation 


3 6 REALIZABILirY CONDITION 01'' THE WIGNER DISTRIBUTION BASED 

WTSaWcPOKS-*" 

For any given time- frequency distribution C(t, W) 

1 f It is possible to find a signal f(t) such that 

f (t) N£^£(t, w ) =C(t, w) (3 99 ) 

then we say that C(t, w) is WD realizable Realizability 
condition for the WD has been stated in property 24 of section 5 5 
Ghd can be derived with the help of basis-vectors in the 
same fashion as in the case of Ambiguity Function [34] « 

Let ^ complete orthohonnhl set of basis funci 

ions Hence a time signal f(t) can be e 3 <pressed as 


f (t) = 


4.0Q 

E 

ias-o* 




(3 100) 
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where the coeffierents are given hy the equation 

4 . 0 a 

=» / f(t) <P^(t) dt (3 101) 



4-00 4-00 

+«» 

y(t-|)e’“^“^ au 

W£^f(t,0) ) 


S aj ' 


„eo l=s„oc 






= S S a.a * 
i 1 ^ ^ 

i (P^(t-l~)(Pj’ (t-|) 

dt| 

(3 102) 


MtOO 


4CQ 

Let ^ dri (3 1o3) 

•^00 


CO) - I E aj* 

£p,,j(t,co,] 


It oan be shown that 


(3 104) 


forms a oomplete orthonromal 


set of basis functions of variables t and W This statement 
can be proved with the help of a previous result namely that 
(t# 0)J forms a complete orthonormal set of basis 
functions [ 34 ] / where 0 ) is the Antoiguity function of 

the basis functions V’j_(t) and (h) defined in equation (3*43) 


Using equation (3 46) vie get 
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^90 J.C0 




1 / 

271 






di.c3(J 


^OO —00 


2n 


/ / D, .(t,w) D „ „(t,OJ) at dw 


UJ 


m^n 


•00 «»co 


4> oo4.<x>f ooL|.oa4,o<4.co 

= ;;//;/ 


iM oou, cicuM^eau.ocu.0* 


, (2ti) 


-j:A^^j(T,0) (C..3) e 


-3 (to- <Jt) 


e'*’^ ^^^“■^'^)do.d(3dTda dtd<t> 


(271) 


^CJO +00 +00 ^00 +CJO +00 

J 1 S 1 J 1 a^^j(t, 0)A<^^^(a,p) -155 

— oo »i«ao **«oo — oo •moo ^co 


^-J (i)(x-u) Q atdW 

+ 00 +O0 +CO +0O 

2^ / / / -p)6(T-a)dadpdcdO 

p^co ^ta *^00 »«oo 


+ eo +»4 

= 5^ (t, a ) (3 ”505) 

forms a complete orthonormal set of basis 


Hence 
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functions of variables t and W This inplies that the time- 

frequency distribution C(t, w) can be written as 


C(t, w) 






(3,106) 


where the coefficient is given by 

(3 107) 

Comparing this with equation (3 104), we find that C(t,t»i) can 
be an auto WD iff 


B 


ij 




for all i,j 


i e iff a » a (3j08) 

til 

where B is a square matrix with (i#J) clement as and a 

til 

is a column matrix with i element as Similarly it can 

be stated that C(t, w) can be a cross V7D of same time signals 
f(t) and g(t) iff 

i “ a (3 109) 


where elements of the matrix b are given by an equation similar 
to equation (3*1o1 ) 

« / g(t) cp^(t) dt 


M^OP 


(3 110) 
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ractorabila.-ty of the square matrix g into the sirrple or ouber 
product of two column matrices §. and b as given by equation 
(3 109 ) irrplies and is implied by unity rank of i [55] 
Moreover by equation (3 107 ) 

4oo 

®Ji 0^(t/ w) Dj^(t,w) dtd(*5 

« ~ fj GV(t, w) (t^w)dtdto (3 nil) 

Using equation (3 16) 

®ij “ ®*i “ c-^Ct, w) 

But C(t, (/)) « C''' (t, w) only if it is an auto WD, Hence we get 
the following result CCt,(i)) is a WD iff rank of B is unity and 
moreover if 

I = (3,112) 

then C ( t, W ) is an auto WD 

3 7 RBAhIZABILITY COMDIllOlSt FO R LINEA R IHTBGRaL TRANS TO RMATIOH 
A linear system in general has following input-* 


output relationship 
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g(t) = / K(t/C) f(T) dT (3 113 ) 

where f(t) is the input# g(t) is the output and K(t#T) is the 
impulse response It can be sho’/m by substitution that the WD 
of the output is given by 

*f€0 4^ 

^ ^ ) dT^do (3 114 ) 

where j^{t#w# t# a) is a two-dimensional WD of the inpulse 
response K(t#t) defined as 

4-co 

^ "In ^ ^ T4^')K*(t-^# T-|') 

^OO ,.«90 

^-j( un-on.) 

In general suppose a kernel X(t# w # T#a ) is given to us and 
we wish to ascertain the invariance of this kernel to the WD* 
The invariance condition will be satisfied if X(t# (i #T/a ) 
corresponds to a two-dimensional WD of some linear system 
inpulse response This in turn implies and is irtplied by the 
fact that the linear transformation of the WD of any signal 
using this kernel is also a WD Hence we wish to find whether 
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^OO ^00 

B(t,w ) 3 ; ; x(t,W/T:,o ) w^£(T,a ) dxdo (3 116) 

mCpO 

correspond to an auto WD of some siqnal or not 

Now according to equation (3*68), B(t/W ) can be an auto WD iff 

+w 

) e do' w g(x)gMy) 

* (3 117) 

foi some signal g(t) 

Now Ii H S of the abo</iS equation can be written as 

*|« C30 oo C0 

^ / X(2^,0),T,cf ) W^^^Ca^a) 

-.00 „oo (5 1l8) 

Let ; dw = \(x#Y^T^a ) 

^00 

^ ^ ^(x,y/x,a ) £('^/o ) dtda 

..M^OP —OP 

Let us eiqjand ^-^£^£C''^/0 ) and X ('C^y/T/ o ) using the complete 
orthonormal set discussed in section 3 6 and 

defined by equation (3 1o5 ) 


(3 119) 
(3. 120) 
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vV£^f(c,a)= s E ^ ^ (3 121) 

i=,*CO j £5,-C0 
Jjp^ ^CO 

and X(x/Y,T, 0 ) = E E (x,y) (t, cJ ) .(3 122) 

3^2>*cd j :=W)o 

where the coeffrci<-nts aj|^ s arc aefined in equation (3 lOli ) 
and coefficients C. , (x^y) are defined as 

1/ j 

Cj^^j(3c,y) // X(x,y,T^0 ) i>i^j(wa) dxdo (3 125) 


TB(x,y) = E S aj^aj- (y,y) 


Thus B(t/ 0 )) Can be an auto WD iff 


J j aj* (x,y) « v(x)v"(y) 


(3.124) 


(3 125) 


for some function v(x) . 

A sufficient condition for this to happen is 


^i,J ~ ^i^^^ 

for some functior 3^^ (x) . 


Sj«(x) 

But the converse 


.(3.126) 

is not proved so fir 




I 






— .r-**** ►.*. -«t# 

A* 91901 ; 
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Special Case 


Let X(t; w ,1,0 ) xCt-i, w -0 ) 


(3 127) 


■Iwo-dimensional WD of the 
which m'=>ans if X(t,w,T,o ) happens to be a /impulse response 

of a linear-system then that lineai system is time-invariant 

Hence equation (3 116) becomes 




B(t, 0)) = / / x(t-i:/^ -a ) ) atdo 


(3 12S) 




Let t-T s= T], 


(1) 


-a 


— w 

A 


B(t, 0)) 


I f xihr V w. ^(t-rj, 0)- w) dhdw 

X, r 

•Oft --0& 


• Tg(x,y) « / 1 


/ x(n, w) u 

•-CO 


« ) 




Now 


Twj f(^'- ^ a» 

^00 -p>o j ^ ^2^ •“ ^ ) 

/ / ff(oi_^) f) e 




«^0O 
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Let = q 

Above integration oacomes 

H-“ 0(^-1) (p-q) 

/ I P(p) F”(q) e ® 


= aJ // 

-CO 

c: e 

Tg(x,y) - 


P(p) e rMq)e dpdq 

f(x-'?7) £Uy-^) 

^00 

/ ; XO?/^) f(3o-r/) f*(y-??) drjd'P 

•-•CO •^oq 


, B(t,W ) is an auto WD coriesponding to some signal iff 


»|«CO 1^0# 

/ / X(r7,^ ) £(x-r3) f’«'(y-Tj)d77d’J' *= v(x)v'My) 

-w -<?• ^ (3 '129) 


for some function v(x) 
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3 ♦ 8 PISCRCTS WD FORfUL 

Com:)utF>t.ion of the V© values involves evaluation 
of a continuous integration which ultimately will be done by 
using numerical methods This inherently involves seme error 
and a lot ol computational efforts To overcome these draw- 
backs we would like to define discrete WD in a manner analog - 
ous to the definition of the continuous WD This new defini- 
tion will use discrete time signal (v/hich can be obtained by 
sampling the continuous time signal) / and all the integrations 
will be replaced by summitionS/ simplifying the computat\ons 
in the process This definition as gi^ren by Claasen and 
Mecklenbrauker [55] for two discrete-time signals f(n) and 
g(H:^is as follows 


1 “> 

u) =22 f{n+k)gMn~k)e~^^^^ 

]C=3^ 00 


(3 130) 


where o is a continuous frequency variable To sircplify the 
computations furthei even the frequency variable should be 
discretized But this should be done in such a way that all 
the results and properties of the continuous WD should remain 
intact In the following analysis the criterion used for 
sampling the continuous frequency variable w is to have 
sampled basis vector j analogous to the continuous 
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basis vector defined in (3 103) This treat- 

ment is also an extention of a similar treatment on the 


Ambiguity Function [54] 

Iiet f(k) he a discrete time signal defined for 
k - -K, -(K-1), 0, — f (K'-l)/ K and 

f (fc) * 0 for jkj^ K 

Let ^ Sampled orthonormal basis set over the interval 

-K to i-k# which implies that 


o (k) = 0 for I kl > K 


(3 131) 


+K 

and S <P^(k) <P* (k) » (3 132) 

k«-K 

Therefore any discrete time-signal f(k) can be expressed as 

+K 

f(k) = S fortk|<K (3 133) 

i— — K. 


To obtain sampled basis vector start 

with the corresponding definition in analogue case given by 
equation (3 1o3 ) 

•f® 


•J90 


56 


Let « 3 


D^j(t,a)) =2 J tp^(3) cPj(2t-.3)e'’J^‘^ci3 (3 154) 


Observing this equation iLet define the sanpl<=‘d basis vector 
Dj . (n,m) as follows 

ly j 






K 

S 'P^(}c)<P* (2n~k) 
k«-K 


2 ?jTmk/^ 

6 

(3 135) 


which is obtained after replacing t and of equation (3 134) 

^ , 27Tm respectively# where P is an integer This means 

by n and ^ ^ 

^ 2t1 

that Sample points in the frequency domain are distanced 

radian J 

Let n and m take values frorm -N to N and -.M to M respectively 
where N and. M aie some integers 

Let us def jne discrete inner product of (n#m) and Dp^g(n#m) 
as 

N M 

A, . „ „ = S Z D, (n#m) D* (n#m) » (3 136) 

ni=^-N m=-M 


Using equation (3 135) and rearranging we get 

N K K 

A , = S S h ‘P. (k)<p!!(2n-k)cp<(h)q> (n~h) 

%J/P/q i J P g 

n.=s-N k«— K h=-K 

M 

„ (k-h)4Triify'p 

h 0 
m5=-M 
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The summation under m can be shown to be equal to [ 54] 
sin { (m+I*) (k-h) 4n/p I 

Let M « K and P = 2 (2k+l) 


Summation under m *= s (2 ki-1)6^ 

rkhiT 

£or I k j f j h j^K 



•I I\ 

H 

Aj j a 

j /P/ Q 

(2k+l) S 

1 V () )tp*(k)(p;(2n^k)^p^(2n-k) 
1 p J q 


ka-K 

n=5-*N 

K 

2N~k 


= (2K-^1) S 

S 

'/'i (k) <P^ (k) <,<j (n’ ) 

k?=-.K 

n’ =s-2N-k 

where 

n* ~ 2n-k 


Now if we put N= 

eK v/e find that for any value of k between -F 

to K, n* always 

covers the range -K to +K 




(3 137 ) 


Thus the basis vectors j oithogonal to each 

other, a result which we were aiming at 
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, with N a K, M - F* p « 2 (2k 1-1) we have 5d . (n^m) C as oornplote 

L / 

orthonormal set of basis vector®* for the disctete WD This means 
we should have {2 k rl) points in time domaan spaced equally and 
(2K^1) points in frequency domain spaced at a distance of 

With this frameworV the results derived in 
sections 3 2 to 3 7 follow for discrete WD also 


3 9 FIST COMPUTaIION OS’ THE WD 


As has been mentioned in the last section for the 

purpose of actual comoutation of the WD on a digital computer 

we have to sample both time and frequency domains Tor the 
purpose of fast computations an ejqpression of the WD similar to 
the Discrete Fourier Transform e^spression can be derived starting 
from the analogue signal The Fast Fourier Transform algorithms 
Can then be used directly 

Consider an analogue sigial h(t) and its Fourier 
Transform H(f) as shown in Pig 3 1(a) To sairple the signal 
in time domain h(t) is multiplied by a pulse train (t) and 

to Sample it in frequency domain the Fourier Transform of the 

resulting signal is multiplied by a frequency domain pulse train 
/N^ f) The time periods of "^0 '^'l 
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respoctivcly. In order to have minimum aliasing ellect we 
mu^t have ^ number of time domain 

sanples which approximate the signal h(t) (i e h(t) is neg- 
ligible after N samples) Pig 3 1(b) to (d) show the sampling 

The final signal 'h(t) can be written as 

h(kTQ) 6(t-.]cTQ-rT^) (3 138) 

It can be argued that the W D of 1i{t) is nonzero only for 
t « ^ k— 0| —— — 

and 0) « 

whore =* 

Also the WD Is periodic in both time snd frequency domains 
with periods and ) respectively Let us try to 

compute the WD value at point.j 

t — k Tq ^ ~ D/1^2^ f (N— 1) 

W = 2Tinfu^ ~ 0,1,2^—, (N-1) 

By looking at the graphical picture and applying the defini- 


op eration graphically 
H"* N-1 

^(t) » S S 

k=0 


tion of the WD we get, 
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Wf ^( t - kTo / « -2nnf^) 


+03 


= S "h(kTQ+pr^) ‘h’ (]a‘Q-pTQ)e-^ 


+00 

» S ‘hCkT^+pr^fh^ <]cTQ-,prQ)e''-^ 

P=3^CO 
I i 

&inco h(t) is periodic we have 

hCkT^) ='h(kT^ + rT^) r = 0,j[,l^+2^ 


Dquition (3 139 ) can be written as 

W^(t » IcIq/ w « 2Txnfj) 

!-«■ N-l 

- y z ^ hCkT ^ aT ^) h ’ 

r«-~ laO 


^ •» N-1 ^ - .14711 11 

« S 2 h(kTQHTQ)h*(kl^-lTQ)e ^ 

ra«.e« 1=0 


-j47Tnr 


N-1 - i47tnl 

- 2 hdcT^HT^) h*(KT-lTQ)e 

laO 


^ 4Ttnr 


(3 139) 


r=— ** 
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The summation under r contributes a constant indep^’inent of k 
and n and hence can be neglected Also 

C'v) 

hCraT^) = hCm^T^) (3 UO) 

where m = o,_Hl,+2, 

and m* ® rn (mod N) 

Wg get 

Wf^,(t m= 2nnf^) 


N-1 

S h (k(3J>l) Tq h^ (k01) T^ g-j^Ttnl/N ^ (3 141) 

1-0 

where indicate mod-N addition and indicate mod-N subtraction. 
Tor the sake of sinplicity in writing we can omit ej^licit 
mention of and and write 

M-1 

W^^(k,n) := S h(k 0 1) h’' (k ^ 

1«0 

Let h(k ^D’^^hCk 01) = g(k,l) (3 142) 

N-1 

,Wj^(k,n) « r (3 143) 

1«0 
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To make this e^qpiession look like Discrete Fourier Transform 
expression let n = m/2 

N -.1 

w^^^(k,m/2) ~ s g(k,l) (i 144) 

1=0 

(3 144) 

yor a const ant k / looks exactly like the Discrete Fourier 
Transform equation Hence to comnute the WD values at points 

t = kl’ / 0) = , m « 0,1,2, — (N-1) 

^ 2 

we can use the Past Fouiier Transform (rPT) algorithms provi- 
ded we compute the value of the function g(k,l) for 1 = 0,1, 

— (N-1) before hand Thus at time t = kT^^ we can get the 
VJD values^ frequencies from zero to 7T(N-l)fj|^ at the Interval 
of Tifj^ radians The WD values for further frequencies can be 
obtained from these WD values by proper extrapolation as can 
be seen from Fig 3 1(d) The same procedure is to be repeated 
to obtain the WD values at other time instants Hence by the 
repeated use of the PFT algorithms we can obtain the WD values 
at all time and frecjuency sanple points of interest 
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3 10 EX^UyiPLBS 

Ji# W.I M !»»■! 

In this section some exirnples have been worked 
out to illustrate the concept of the WD 

1) Let f(t) = 1 I 1 1 < T 

= 0 I 1 1 )> T 

Then 

0) ) = sinf 2w(T-t)] |t|<;T 

- 0 I tj > T 

This WD has a sin x/x shape with respect to frequency just 
like the Fourier Transform of f(t) The width of the main 
lobe however depends on the value of t Also the is 
negative in certain regions of the t- w plane • 


2) Let f(t) w e 

1 1 1 < 

T 


« 0 

M> 

T 


Then 0 ) ) t= 

sin f 2(w - 
*■ 

0 

w^) (T (t )] 

1 tl< T 

= 0 

|t|> T 



ja> t 

3) Let f (t) M Ae ° 

for all t 



Then W) = |a|^ 





\ 
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This result i^ rruito satisfactory be ORUge this show^’ for a 
stationary signal thf' VflD is independent of t an i is confined 
to the* line O) = o)^ as exoceted 

t 

4) Let f(t) « A^e ^ 

J0> t 

g(t) A2e 

, )t n wo> 

Then w) = A 2 e ^ 2tt6(w ~ - 3 ^ 

5) Let f(t) « A cos ( H 

Then w) ^ jA|^ 

2o((«') cos £2(o'^t + (P)}] 

Thus apart from stationary contributions at frequency + ,the 
VtD also contains a nons tationary term oscillating with 
frequency 2^^ 

i t^ 

6) Let £(t) w Ae'^ ^ 

Then ~ 

This again is an extremely satisfying result because it clearly 
shoves the concentration of the V?D at frequoncy w « 2ci,t which 
is the inotaiitaneous frequency of f(t) 


66 


7 This example illustrates ho\r the VID ch-r^nqes \fxth a Ixncar 
iransform'ition o£ the t-w plane For this consider a 
recLangalar pulse of width 2T = 4 anu its 'JD as mentioned 
in exanrole 1 

£(t) - 1 |tl/ 2 

- 0 pl> 2 

ana £(t, w) = ^ 

“0 | t |>2 

Piy J 1 sho 73 m artplitude contour of the WD for a value of 
1 2 which brings into sight only firsL side lobe Suppose 
we wish to transform the plane so that the side lobe is off the 
w --axis To this end consiuer following Lransformiticn 





— 



t' 


1 

8/7 


t 



1/4 

9/7 


(0 

J 


— 




The dotorminant of the transformation matrix is 1 as requiied 
for the invariance condition The t-axis and w -axis are 
transformed to straight lines 9t‘ = 8W and t‘ » respec- 

tively. Fig 3 2 shows the new anplitude contour for the 




i 
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same amplitude value The nev/ signal h(t)/ having the trans- 
formed WD can be obtained by using the results of section 3 4 
and [ 30] Here 


a 




4 


7 


b = 


2(det D) 



k « dot D = 1 


and 


or 


f^(t) » f<t) 

2 

g{t) » f(t) 

h(t) = r T “^ { G( « ) i 

H«o) = f(t) at 


-j 


i8W 

7 


w 2 ti e 




7^ 


s/2n 


+ j S 




j 8 


( 2^)5 

\/2n 


where C(x) and S(x) are Presnel integrals* 



CHrtPTSR 4 


APPLICATIONS OP PHD WIGNER DISTRIBUTION TO OPTICS 

4 1 IWTROl^GTOT 

One of the reasons why the WD has attracted 
atbcnLion of scxcntists workang in the field of signal pro- 
cessing is its successful u^’e in optics PapouLis introduced 
the use of Ambiguity Function in optics in 1974 [36] 

Because of the similar ity m the definitions of the Ambiguity 
Function and the tID the logical step was to extend the 
analysis, with the help of WD This work was done by Bastians 
[3] I [ 7 ] Bartelt and Brenner [5], [6] Brared/ Gram 

and Schenzie [4] and others They also introduced many more 
new applications* This chapter briefly describes some of 
the application areas to bring out the usefulness of the WD 
in optic'^’ and to give an idea regarding the use of WD in 
signal processing 

4 , 2 TERMINOLOGY 

Before presenting the applications to optics it is 
essential to introduce the terminology used In Pourtor OTpdiios 
An optical signal can be described in the space domain by its 
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complex arnplitude f(r), where r denotes a vecter with space 
coordinates (x,y) for a two-dimensional signal Equivalently 
in frequency domain we have the Fourier ‘transform of the 
complex amplitude defined as 


r (m) » / £(r) ej<p {•-j(o^r)} dr 


(4 1) 


where w is a t\;o -dimensional frequency vecter/ / dr stands 
for IS dxdy (for a two-dimensional signal) and u r indicate 

«MtOO 

the dot-product of vectors and ^ bimil arly a ] inear 
system is defined by its impulse response g(r/P_) or its 
Fourier Transform 0(^/2.) defined as 


G(W/a) 


it? 


dr d^ 


(4 2) 


The outpuit fQ(£) of such a linear system with input f^(r) 
is given as usual by 


fQ{r) =f / g(r/P) f^(P)dP (4 3) 

The WD defined with this terminology reads 

W^^fCr/W) ^ s f(r h f"(£ - I* (4 4) 


For a linear system defined by equation (4 3) the WD of the 
output is given by 
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Mf ,£ (£<e.) 

o' o 


g,g ^ 


where Wp i'^fO) is the WD of the input signal and w 
■^i 

is the two-dimensional WD of the impulse response of the system 
defined as 


Wg^g(r,oi.,a/2.) * g(r »• f~ / E I* 

if 3r^ P ^ 

»g a " i" ^ - j (^1 t;'- a dr'dP* 


. (4 6 ) 


The function ^ considered as response 

g / g 

in the space-frequency domain to the input signal 

Of course this signal does not have any physical origin but 
only mathematical convenience* 

In the next secticn we consider a few elementary 
linear optical systems* In each case the WD of the impulse 
response of the sy stem is given along with the VJD of the 
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output, signal It will be seen that relating the WDs of the 
input and output signals is much more convenient than relating 
the signals directly either in space domain or in frequency 
domain 

4,3 EXAMPLES 

1) SPREADLDSS SYblEM LENS 

Let the spreadless system be represented by 

Wg^y(rWi>.,a/2.) = 6 (r-o) (4 8 ) 

^ 0 ^ 0 47X i'^1 

. (4 9) 

Equation (4 9 ) shows convolution for frequency variable ^ and 
multiplication for the space variable r# a r suit quite 
convincing for a spreadless system* In the special case of a 
lens with a focal distance f we have 

m(r) » 0>q> { 


(4 10) 
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Here the wave number k it. defined as 


k = 27 iA 


whore /v is the wavelength of light 


= 6 6(w - o) 


and hence 


Wf f (r,w) = w. . (r/y-| r) 


2) SHIFl-INV/allAl^r SYfalEM PREZ: SPACE 

Shift-invariant system is the converse of a spreadlcss 
Let 

represent the shift- invariant system 


and Ww ^ iRtH) =* Wv, ^ -p (P»<^')dP 

Cq/Xq n,a ^i>^i 


(4 11) 

(4 12) 

(4 13) 

system 

(4 14) 

(4,15) 

(4 16) 


This shows a convolution for the space variable and a multipli 
cation for the frequency variable a result as per the expec- 
tations In the special case of frv^e space over a distance z 
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we hstve 

H(w) « C3^ {(-j |w| 2)} (4 17) 

where k is defined in equation (4 10) 

~ <5 Cr - ~ w • p;6(w-*o) (4 10) 

and hence 


W 




(4 19) 


3) FOU^ILR TRANS iORMLR 


Por a fourier transformer whose ii-npiJilse response 


reads as 


9(r,£) « ^ exp{ (-J/3r £) } 


(4 20) 


where j3 is a constant# we have 


(j) 


Wg^g(r/O#£,o) * 6 0r - o)6(~ H £) 


(4 21) 


and hence 


-W 


Wf (r/W) * 


(4 22) 
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It can bo seen that the space and frequency domains are 
interchanged as should be for a Fourier Transformer 

4) MAGNIFIER 

For a magnifier the impulse re so on se is 

g(E/£) « t 6(tr P) (4 23) 

where t is a constant 

* « 6 (t.r-.£)6(^ ~ 0 ) (4 24) 

and hence 

f (r/S) « W . . (r, tr, J) (4 25) 

5 ) FRE SND L ‘ S DI FFJIACTION 

For a Fresnel.* s diffraction we have 

(Jkz) exp * (4 26) 

where z is the distance and k is defined in equation (4 11) 


Wg^(^(r,W/P#2.) = (^Cr^) - « ) 6(|(r-P) -o) ,(4 2?) 

' X a 


Q/g 
and hence 
W 


£ ,f (£-1) » C -t') 

o' o i i 


(4.28) 



CHAPTER 5 


A RECEPTION SCHh,MD USING THE WIGEWER DISTRIBUTION 

5 1 INTRODUCTION 


In this chapter a new reception scheme based on 
the WD is proposed. We consider a transmitter ti •'nsmitting 
a series o£ pulses which reach the receiver after a series 
of reflections from a target as in the case of rsdaar signals 
or as in seismic applications Thioughout this discussion 
we assume that the channel is distortionless and the noise 
IS additive white gaussian Thus the received signal r(t) 
can be written as 

r(t) « S a f(t«»t.) + n(t) (5 1) 

1 ^ 

where a^^s are the reflection coefticaents# x^s are the delays 
involved/ f(t) is the transmitted pulse and n(t) is additive 
white gaussian noise The receiver is required to find out 
the reflection coefiicients a^^s and the delays x^s Since the 
noise is white gaussian the correlator (or the matched filter) 
reception is the optimum linear recexption scheme as far as 
the signal to noise ratio (SNR) at the output is concerned 
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But this reception scheme gives rxso to the sxde-looe problems 
tind may shield certain return if the reflection coexficient 
associati d with it is small The analysis included in this 
chapter tries to investigate whether this problem can be 
overcome by using a reception scheme based on the MO 

6 2 l)^iJVAxION BDIIINP USING THE SIGNER DISTRIxBUION IN 

Throughout this analysis we assime that the pulse 
f(b) transmitted by the transmitter is a pulsed cnirp (or 
linear FM) signal with carrier frequency , i e, 

Jw t 

f(t) « A e ° for 0 < t < T^ (5 2) 

w 0 otherwise 

where is the duraticn of the pulse and A and a are constants 
We assume that >> 2 a which is the approximate 

bandwidth of the signal f(t) 

f(t) bo apnaidorod ae a product of two signals, 

g(t) = e e ° for -«iK (5 3) 

and h(t) =» A for 0 <; t < ^'q 

ss 0 otherwise (5 4) 
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Let 




f(t) 



4- 

— > 

g(t) 

and 

%h<« 


-> 

h(t) 


It can be shovm that 


Wg g(t, W ) = 2 tc6(W~o^ ^ 2at) « (5 5) 

Using equation (3 61) we get 

Thus the WD oi the signal f(t) is same as that of h(t) except 
for a shift of the -axis at each time instant t It can be 
shown that 

= 4At for 0 < t <1^2 

sin(2W (T -t)) 

t= 4A(T^-t) T“ ^ 

= 0 otherwise (5 7) 

A pictorial view of Wj^^j^(t/ ^ ) at some time instant twt^ is 
shown in fig, 5 1(a), The WD of the signal f(t) will be 
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©vcactly like this except i-or the fact that at any time instant 
t “ t^/ the entire distribution is shifted by i-2ut^) along 
the positive (»' -axis as shown in fig 5 1 (b) Thus the ♦cntral 
lobe of the WD of f(t) at time t will be at a frequency 

*** ^Iso adding the distribution values of f(t) 

along line (i> l 2u,t + b (where b is a constant), is equi- 

valent to addiny the distribution values of h(t) along a hori- 
zental line w a b The WD of h(t) reveals thit at any fixed 
time instant t ( < T^/2) the WD has central lobe at o and 
the first side lobe at ^ (approximately) This moans 

that if we add up the distribution values of h(t) along 
different horizental lines b, we ejqpect a sharp peak at 
b = O and since the secondary side-lobes do not f-^ll on same 
horizental line, we expect that as b increase (or decreases) 
the summation should decrease contJLnuously witha-it any prominent 
side-lobes. This, if it works, will be in contrast to the 
correlator (or matched filter) output which shows prominent 
side-lobes and hence can be used in a reception scheme With 
this motivation the following sections try to work out details 
of such a reception scheme and conpare its perfoimance with 
that of correlator leceiver 

5 3 RSCEPTIQN SCHSME USI NG TI-I B CORI ^L ATOR 

Let US for the moment assume that the received 
signal r(t) given by equation (5 1) consist of only one return with 
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unity reflection coefficient and time*-delay t Also let the 
constant a of the equation (5 2) be such tliat the signal f(t) 
will have energy E 


A - 


Thus the received signal r(t) becomes 


(5 8 ) 


r(t) = gJ'- I- n(t) for i < t (ciT^) 

o 


n(t) oth n/i^e 


(5 9 ) 


Our objective is to determine the time~delay t binca the 

noia. is additive whibe gaussian^ the correlator (or matched 

filter) recejition is the b -n^t as foi. as bNR at the output is 

concerned In this scheme we correlate the received signal 

with the reference signal f. (t) = f(t-.t ) and this correlation 

o 

is carii-^d out at all time instanus to Thus the output C(t) 


of the correlitor is given by 

t h r 

C(t) = 


O 7^ 

/ r(r)) f (g-t) dn 

t 


(5 10 ) 


fhe outout C(t) consi'jts of tv/o parts 


t IT. 


Signal part C„(t) = S f(7])f (ri^t)m for (*1-1^) < t < (x iT^) 
3 V O ^ O 


a 0 otlierwiao 


(5 11 ) 
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t -r T 

and Noise Part ^ n('iT)f (T)~t) d?! 

t 


(5 12 ) 


The signal part C (t) is maximum when t ~ u and is equal to 

s 


G^(t) 


max 


G (t - t) => E 
s 


(b 13) 


By observing the outouh C(t) v/e wish co find out the value of 
X Since the signal part of the output O (t) is maximum and 
real at t = T, we musL. concentre. te only on the real part of 
C_(t) Redefining C^(t) as the r^al part of the signal part 
Oi- the output, wc get 

t '-T 

Cg(t) R P S f^T7)f (7?-.t)dn for (c-T^) It (T+T^) 

t 

- 0 otherwise (5 14) 


where R P stands for 'Real pairt of » Since we are mainly 
interested in studying th.^ side'^lobc behaviour of this recepition 
scheme wc can assuima 

T = 0 (5 15) 

By carrying out the intogration of equation (5 14) it can be 
shown that 



S4 


C (t.) = sin{u.‘t(T -t)} cost o) t + C\-tr } 

^ ^ O O 

for 0 

~ rrTrV^ sm { U-t(T rt)} COs( t + C.tT ) 

X wuTI. o o o 

o 

for --i < t < 0 
o — ~ 

= 0 ofnerwise (5 16) 

The derivation of thio result is given in Ippendis i 
The envelope A^(t) of the output Cg(t) is given by 

AQ(t) « I sin tat (T^-t) }[ for 0 <t <T^ 

“ T“tt^|t| I ^ ’•‘^o ^ ° 

= 0 otherwise (5 17) 

The noise part C^^(t) of the output is a stochastic process 
with gaus&ian distribution The variance a of this noise 
can be shown to be 



(5 18) 
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wnero N/2 is the noise spectral density of the input noise 
n{t) (ser <i!^pendix \) 

&NR at the output when the signal part of the outout is 
xnaximum is given by 

(bNR) , . = w" (5 19) 

^ correlator N 

staxted this section with the assyrrption of only 
one r-turn with unity reflocticn coeflicient For a received 
signal consisting of many returns with different rcflecbion 
coefficients as given ty equation (5 1), the correlator output 
v/ill be the superinposition of the outputs of the individual 
returns The outputs of the individual returns will have 
peaJes corresponding to the delays of the returns and the 
airplitudes of the peal s \';ill be proportional to the corresponding 
rofioction coofficients Ho\/ever each of those outputs will 
have secondary side-lobu-s apart from the main lobe llcnce the 
superirrposition of the individual cutouts will still have peaks 
coriesponding to each return but thu time instants of the 
peaks may not exactly match the delays of the. returns 
Similarly the amplitudes of the peaks may nob be exactly pro- 
portional to the reflection coefficients of the corresponding 
roturnis » 
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In sectxon 5 6 simulation results of the correlator 
output havo been presented assuming two returns with variable 
delays and reflection coefficients 

5 4 RECEP TION SGHGIIE_ USIdG^ PIS IRIBUflCW 

The reception scheme using the WD consists of taking 

the cross WD of the r-cciv d signal r(t) and tho reference 

signal f+. (t) ?= f(t-t ) and then summing uo the distribution 

values along the lin^ o' + 2a(t-t^) of thi t~6) plane This 

procedure is repeated for all time instants t^^ and the instants 

proVDd^ a measure of the 

and amplitudos corresionding to the p^aks in th-^ outputyref~ 
lection coefficients and the time delays as in section 5 3 
we assume the received signal to oa consisting of on^ r-^tum 
only with unity reflection coefficient and tima-delay v At 
the output of the receiver ’i/e expect a peak wh^n which 

should fall rapidly as t^ increases further 

The output W(t) of tho WD based receiver is given by 
t+r^ 

W(t) = ; ^ + 2Cu(0-t) ]d./ (5 20) 

t 

when f^(3c) = f (x-t) . 
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Using -actuation (5*9) wo get 
t+r 

O 

w(b) =- / W ^2c(0~t)]dr] 

t > * t 

for < (x+2To) 

t KT 

o 

” f *■ J otherv;isw 

^ (5 21) 


This can be written as 


t+T. 


t-lT 


W(t) = 


/ =o> ^+2a(7/«t) ] d7? + ^ \ 

t V t ^ ' t 


o)^i-2a(r)^t)] dr) 


for (c-2To):^’t <('^‘*’ ^T©) 


t i-i 

= / W = 0) +2a (77-t) ] dU otherwise 

^ * (5 22) 

The output v?(t) consist of two parts 

ti-T„ 

o 

Signal Part W„(t) « / W. . = w H 2u(?-j-t) ]dr) 

S Jw X o 

t ^ ^ 


for (T~2To)<t < (tt 2 Tjj) 


w O otherwise 


(5 ?3) 



t+1 

o 

I 

t 


^8 


and Noise Part ® 


W ^ C 0 ) = W h2a(r)-t) ] d?7 
riy o 

(5 24) 


The sigiial part Wg(t} is m'^ximum when t = t and is equal to 


W3(t) 


max 


W„(t-T) = ET 
s o 


(b 25) 


The derivation of this result will be presented snortly 

Again since the maximum value.-* of W^(t) in which wc are interested/ 

3 

is real/ wo should concentrate only on the real part o£ /L(t) 

s 

Redefining Wg(t) as the real part of the signal part of the 
output/ we get 


t+T^ 

o 

Wg(t) »Rp / W^^[f)/W » o)^-i-2a(r}-t) J dTi 

t 

for (T:"2a;o)< t < (T-rat^) 


= 0 otherwise 


(5 26) 


As is done in Section 5 2 let us assume T « 0 The value of 

W_(t) at t = t can bo obtained by first computing the WD/ 
s o 

Wf ^ (t/W )» Referring to Tig 5 2 and using the definition 
oi tliG WD it can be seen that 

^ (t/ o) - 0 for t<t^/2 or t > H t^/2 

* t 

o 


(5.27) 
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For ■b^/2 < t < J? l- t /2 

O — *”0 o 

Wf f (t,u ) = ! f(t4) f an 

' _<» ° 

4-#ci 

= ; £(t}^) f*(t - I - t^) an 

Let n ss S -* t 

o 

4- 00 

' ^ "t 

(t,w ) » s f(t - ^^ + |) f*(t ~ -2^- - §• ) e"*J“ (G-t^) 

j(i> t t 

= e ° ° (t - -gP- ,w ) for t^/2 ty2 

(5 28) 

Hoiice summation of the WD values £ (t, <^ ) along an a^cis 

^ t 

^ ~ o ' is equivalent to the oumniation of the WD values 

W- f.(t/ O' ) along an axis O' -o' j-2iju(t-t /2) with t going from t /2 to 
XfX. 00 O 

n?o+tQ/2 and with eaoh WD value being multiplied by e ^ As has 
been e^tplained in section 5 2, this in turn is equivalent to 

the summation of tho WD values ^ ' defined in equaticnj 

. . ' with . ”^0 

(5 7)/ along a horizontal linow- -Cut^t going from t^/2 to 

lO^ * 

and v/ith each WD value being multiplied by e , where 
0) * ® Wq + 2tt(t-t^/2) 
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Hence 


Wg(t) 


‘ Vh C='' J' 

t/2 


j[0)^ + 2a(x-t/2)]t 


(5 29) 


By using equation (5 7) and carrying out the above integration 


\/e get 


Wg ( t) =0 for t ^ - 2T^ 


E(2T+t) 2 r 2 

J sin(?atTo-tt ) + — sin(OutT^) sin(ab r^-Kit ) } 


2T t~t 
o 


cos (^^t) 
o 


E ( 2T +1;^) 

-Vn-r'S' <2os (2otT -a,t^)_ — cosCcutT )sin(ctT -at^)} 

^ slnCw^t) 

for ~2T^ < t < -*T 

o — o 


cos (u'tf^-at^) s±n(tfta^) t sin (2cutT^^at^) 


^ 2 2 
2TQa t 


sin(ct'^) sin(2tctT )} cos t> 


+ (2a'tT -.ut^) - 2 sin(at^) cos(2atT ) 

■^o'^ o 2 ^ u.^t 

o 

- sin (ci-tT -at^) sin(ati )} sin (^^t) 

UJUu O O O 

for - T < t < Q 

O — * 
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s(Vt) ^ 

{ 2*ruir"“' (2o,tT^„at )} cos (W^t) 


cos ( 2o,tT^-.o.t^ ) w — 


'2j^ 

^TJX t 

o 


2 ® 2 
sxn at ~ ^ 

o 


sxn( ^'^t) 


for 0 < t < f 


2 F ? 

{ V TT — Sin(2ata' -at ) — • %— sin{2atT -at ) 

^o ° 2T a^t" 

o 


sin (at )} cos( 


E(T -t> 

+{ — cos (2atT -ct ) — sln(2atT -at ) 

° 2T a^t ® 

o 


GOs(ct )) sin (^Q^) 


for T < t < 2T 


^ Q for t > 2T^ (5 30) 

Derivation of the above results is given xn Appendix B It is 

interesting to note thaU W (t) is not symmetric about the W (t) 

s s 
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axis The envelope q£ the outout can be oo tamed by squaring 

the coefficients of gos( and sin(t'^t) terms adding them 

and taking the square root The value of /7_(t) at t =* 0 a.s 

s 

maximum and equals EJ?^ as mentioned in equation (5 25), 

The noise part of the output Wj^(t) is a stochastic 
proceso with yaussian distribution The variance of this noise 
can be obtained by applying the definition of the variance to 
equation (5 25) Rewribting equation (5 25) we get 

Wj^(t) = ^ “ ^o ^ 

t 

Let X =5 T}-.t 

^o 

Wj^(t) =» ; + 2ax] dx 

0 

~ ; n(tTX f ^ ( ±~^) dhdx , (5 31) 

X~0 ?7=-i0o 

T^ ( oo i-co 

W*(t 2 )]=S f / S f n(t^+x-j^)n'(t 2 ly 13/2) 
X=0 y=0 


f U^n/2) f(y^p/2),e 




JV 


J2ay|3 > 

i -» 'qJ dj3dndydx 
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Sinoe nCt) is wnite noxse wit-h spectial density w/2 we havs 

E [n(tj_) J= I (5 32) 

1 T <-oa Hoo 
O O 

L CWj;j(t3^) Wj^(t2)]= J J f J I 6 (t^+X+|-t2-y-^) 

X-O Y «=^> ))=:>« oP=!-«i 

f (x-n/2)f(y~p/2) dpdtidydx 

T i 
o o 

=• ; J 

x=o y=o 


+00 

/ 


n 




/7 Qt^OO 


^-J2ax^ ^J2ay(2ti4.2xM)-2t2-2y) 


Wow for given x,y,t.^/i 2 e cist values of T} such that 

(x~ 2 ) Q’hol (2y+t2"t^-*x-J7/2) both Ixe in the interval [ 0 , 

Lot R be that range of 77 Hence whenever rj belongs to R \/e 
can substitute the egression of thi'=* function f given by 
equation (5 2) 

. E [w^(t^) w (t 2 )J “ ; ^ T 

X«0 7(CR ^ 

jw (?y+t «t.-7)/2 -x) jOi ( 2 t^l- 2 x- 2 t 2 - 2 y) jaCay+t^-tj^-'Z/a-x) ^ 


Q~J SCbXT ? 


j2c y( 2 t^+ 2 xr'n~'^t 2 - 2 /) 

© 


d7/dydx 
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0 Bn 

j a[ (t 2 -t^ ) ^- 2 x(t 2 -t^) ) ] 


= ; / p~ 

0 

e 


x-o y=o 

(t -t ) 



e ^ •" d//dydx 

(5 33 ) 


Thus E [ Wj^(t^)Wj^(-b 2 ) ] IS a function of (t 2 -t^) only Hence 


the noise is stationary as expected fo commute tho 

variance of t^ =2 t^ in aquation (5 33) 

T T 
0 o 

^ E V^(t)] = / / ^ dr/dydic (5 34) 

x~o y o OLR ^ 


The region R can be found by separating the intervals of x and 

y into smaller intervals and then carrying out the triple 

intejration of equation (5 34) This can oc done and the 
2 

variance 0 of the noise can be shown to be equal to 


0 


2 


bENf^ 


(5 35) 


Derivation of this result is given in Appendix C 


The Sjsir for maximum output 

(W (t) )^ 

e rnaxt__ 

0 ^ 



(SNR) 


(5 36) 


WD 


/-T 2„,2 

6L 


SENP 


6 E 
5 


Comparing this valu^ with the SNR value ox the correl- 
ator receiver (equation (5 l9))/ we se'=‘ that there is a de- 
gradation by a factor of 0 6j\ degradation of -2 2184 dB as far 
as faHR performance of the WD based recxolvar is concerned 

For received signal consisting of more than one 
returns the output will be the superiiTiposition of the cutouts 
of the Individual returns As in the ease of the correlator 
receiver output, this superimpositian also will have pealcs 
corresponding to each return but the delays and the amplitudes 
of these pealcs may not be an exact measure of the delays and 
the reflection coefficients of the returns Section 5 6 
presents some of the simulation results assuming two returns 
v^ith Variable delays and reflection coelficients. 


5 5 WD BASED RECCIVER AS AN EQUIVALENf LINEAR FIIiTER 

Prom equation (5,20), which gives the outnut of tho 
WD based receiver, it is clear that this reception is linear 
and time -invariant in nature Hence it can bo replaced by an 
equivalent linear time-invariant filter 



97 


The impulse response o£ such a filter can be 
obtained by considering 


r(t) = 6(t) 


(5 37) 


in equation (b 20) IL can be shown that the required impulse 
response h^(t) is givon by 


h (t) ^ 0 
w 


V “ ( 2 ^ 0 ^) e ° e 


-jut^ 

6 e 


-s/® ®0 


2T > t 
o 


-.2J? 


o 


-J? 


«T < t <0 

Q 


0 lO)^ 


= 0 




(5 38) 


Derivation of this result is included in ^pendiwC D 

The irq‘)ulse response of the correlator receiver h^(t)j^ is 
given by 


h^(t) « 0 



at 0 


-T, 


o 





<t< 0 




(5 39) 
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Compdrinnj the two impulse responses h (t) and. h (t) it is 

"Vf c 

seen Lhit in the interval ^t <0, both have exactl/ 

same Value except for a scaling factor of and sign 
change of the carrier frequency term Outside this inteirval 
hj^(t) is zoro v/hereas h^(t) is still nonzero on both sides 
upto a length of 

The output of the WD based receiver nov/ can be written as 

V7(t) » r(t)* (5 40) 

As has been mentioned before for the case of single 
return with zero delay and unity reflection coefficient, the 
peak of the output occurs at t = 0 in both the receivers 
In the correlator receiver both the signal part and the noise 
part of the output get affected by the entire inpulse response 
h (t) f However in the based receiver the signal part of 
the output is due to only a portion of the inpulse res].)Onse, 
viz iiBtJulse response in the range < t <0, wh'^'reas the noise 
is affected by the entiie Impulse response h^(t) This 
eJtplains the degiadation of SNR at the output of the viD based 
receiver conpared to tho correlator receiver as mentioned in 
section 5,4 At the same time it is lor^cisoly this ‘extra* 
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portion o£ the impul c rosoonse h^(t), which is ejtpect^d to 
suppress the sidelobes in the output of the WD based receiver 

It is now clear that the V)D based receiver; although 
defined by equation (5 20) , need not be inplemented directly 
as per this equation because this involves computation of a 
large numbci of WD values which is inher ntly i coircplesc and 
slow process Instead implementation using the inpulse response 
h^(t); is more convenient from a practical point of view 

5 6 SIDE^LOBgS OP TliE CORRSLAfOR RECEIVER AND TI-E VTO 
B A5ED l^GSI VLR 

In this section simulation results of the outputs 
of the correlator receiver and the based receiver have bean 
presented Since we aie interested mainly in the side-lobe 
charac tons tics of the oub-jut; the noise has been omitted from 
the received signal r(t) 

For a single return with unity reflection coefficient 
the envelope of the output of the correlator receiver is 
given by equation (5 17) Tho envelope of the WD based receiver 
can be obtained from equation (5 30) Both these envelopes 
have bean plotted for different valui^s of the parameters and 
The plots have been presented in Append! c E 
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It can bg seen from the plots with fi cecl “ 2 Secs 
and variable u that for email values of a tne widths of the 
mainlobe and the first sidelobe are small in the output of the 
WD based receiver The maximum arcplxtudes of the sidelobes 
are hov/ever more or loss sam in both the outputs and hence 
there is no suppressroi^ of the sidelobes in the Wt) based 
roceivor Moreover agti* increases even thxs marginal xmprove- 
ment in the lobc~width reduction vanishes and both tho outputs 
hav(. sxmilar Sidelobo ch ^ract ristics 

A better way to comp tre the tvro reception senomes 
for the sidelobe performance is to actually consider two or 
more returns with different reflection coefficients ior the 
purpose of simulation, two returns have been assumed The 
fir^t r*" turn has zero delay and unity reflection coefficient 
The delay T2 reflection coefficient 32# of the socond 

robum are variables The plots of the envelopes of tho 
outputs for different values of aj and Iwe b^en given 
in Appendix y* 

It can bo seen from these plots that for 32 =”^ ' 

both the outputs have srmilar nature for all values of T2 
Both tho receivers are unablo to defect the second return 
clearly for small 13 detect it with same clarity for 
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largo Tj, For smaller values o£ a2(ap«0,l,0 6), tho perfor- 
mancv^ however dilters Here again when is smaller thai 
tho width o-t- the mainlobe/ both the r'^ceivors fail to aetect 
the second return* But if T 2 is such that the second return falls 
in tho first sidelobe of the outout of tho first return (C2“40^50)/ 
the WD based receiver has better clarity or the discriminating 
ability/ i c tho slope of the rise ana the magnitude of the 
ri o oi the peak in the outout corresponding to tho second 
return arc both lai^o in tho WD based rcceiv^^r When tho 
secozid return falls in the second sidelobe of the output 
(Xg « 70/00)/ simil cO. twsult is oba( rvod In this case the 
discriminating ability of the \JD bas'^d receiver is more pro- 
nounced but at the sam time there is a presence of a small 
spurious lobj in between the two lobes (corresponding to the 
two rc turns) / wnich is aosont in the correlator r^ceiv^r output 

Whan a 2 is further reduced to 0*2, similar results 
are observed* When the second return lies in the second sidelobe 
of the output of the first return/ the correlator receiver is 
almost unable to detect it W based receiver/ however can 
detect it clc arly although the anplitudo of this second peak 
is more tlian proportional to Sj* 
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Overall it IS sO'-n that both tho raception schomus 
haVG similar performance for determining the values of the 
time delays In the determination of the reflection coeffi- 
cient values also, both perform equally v^^ll (or bad) , in an 
ov^-rall sense, ivQ* in certain cases based receiver gives 
better results and in others the correlator receiver 
Howevui as Iftr as the detection of the presence of the second 
return with smaller rt flection coefficient is concerned, 
the t© based receiver has distinctly bettor oerfornDneo This 
IS particularly true when the second return Xiei nil the 
first or the second sidelobe of the output of the first 


1C turn. 



CHAPTER 6 


COHGLUblONS 


This chax:>ter xs aimed at reviewing the significant 
resvxlts obtained during Uie course of this work ana making a 
few suggestions foi the future line of worlv in. the applica- 
tion areas of tho WD 

A couhined time-frequency representation of signals 
and qj^stcms is an extremely useful tool in signal processing 
This is particularly so for nonstationary signals whoso frequency 
content changes with time fhe exanple given in Chapter 1 
clearly indie ^tes the shortcomings of the conventional repre- 
sent jtio ns and brings forth the usefulness of time-fr ouency 
representations A general class of such representations was 
pioaosed by Cohen and has been discussed in Chapter 2 HD is 
a member of this class with properti-^vi more suitable for 
an}.:)lications in signal processing 

Definition of the WD is similar to that of the 
Ambiguity runction aiid hence most of the la ter ‘s ground work 
Can straightaway be extended to the WD Chapter 3 gives an 
exhaustive list of tho mathematical properties of the vTD 
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Most- (but not ill) oI the results look similar to those of 
the ambiguity Function but details diflei \ method for 
fast corrqputation of the WD his also been developed in this 
chapter 


vfD also finds its usefulness in optical signal 
processing Lxanpl s given in Ch''ptei 4 are only a small 
gliinpsa of the whole apolication area currently under active 
consideration Motnods for optical generation of the v'JD have 
been suggestuQ [5] , [6] Input output relationships are 
similified in the WD domain In many of the commonly used 
optical systems like lenses^ magnifiers/ fourier transformers 
etc. Study of the diffraction patterns in. a cortplex environ- 
ment will also be &lmplifi..,d in the WD domain 

\ reception scheme based on the ^VD is proposed 
in Chapter 5 to overcome the problem of side-lone interference 
Details o-f this noM schema for the reception of pulsed chirp 
signal have been worked out Compared eo the correlapor 
receiver this receiver has a degradation of -2 2184dB as far 
as the bNR at the output is concerned The sido-lobe per- 
formance howevoi shows some Inpiovement though not as much 
as was thought of For the computation of time— delays and 
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refl-cction coefficients there is no marked iroprovemont in the 
new loGciver but from the point of view of detection of weaker 
returns falling in the side-lobes of stronger returns it 
clearly shows advantage over the correlator receiver The peak 
in the output corresoonding to such a weak return is more 
sharo in the WD based recaiver ijodicating its better discrimi- 
nating ability There are cas^-s when the correlator rcceivoi 
aomplatrly fails to detect a weaker return whereas the V?D 
based receiver does it comfortably 

Some recent developments in the field of WD^ which 
have not bem touched so far in this thesis / deserve a mention 
at this stage, /iD has been used successfully in the design 
of loudspeakers by Janse and Kaizor [15] opening a vast now area 
of applications in acoustics A systematic method of the 
synthesis of a signal from its v® has been given by Kai-Bar 
Yu [38 ] CLaasen and MocMenboker htvo discussed in detail 
aliasing in the discrete WD [39] A begining has been made by 
Kay c.nd Boudrcaux-Bartels in using the WD for detection 
problem [40] , Approach used by them is quite similar to the 

one discussed in Chapter 5 
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njnjjs J3C0JPE oi! mw 

Any analysis involving nonstation arity merits an 
attontion irom the WD Speech analysis is one such exairiple 
The spcctrogiams used e-<cten.aiveiy in speech analysis are based 
on the assumption ol stationarity over the window length. Such 
an artificial us umption v/ill not be required to be made while 
using the tID Moreover the time and frequency resolutions in 
the spectrogram should both be better with the use of WD 
However WD based spectrograms v/ill need huge amount of conpu- 
tations and fast computational methods must be developed 
One such method has been discussed in section 3 7 but it 
needs to be improved further 

The WD based receivei discussed in Chapter 5 assumes 
only time delay in the returns but no doppler shift An 
extention oi the same reception scheme can be made to accomodate 
this also Such a modification will involve summing up of the 
distribution values along lines with different slopes in the 
twO) plane and looking for the peaks It v;ill be interesting 
to invsvjtigate the SNR and the side-lobe charictcri^tics of 


such a receiver 



APPENDIX A 


A. 1 OUTPUJ? OP IHE CO^’ICjliATOR RECEIVER 

BOR A SINGiX RDIDRN ^ ^ 


According to equation (5 14) the signal part of 
the output of the correlator receiver for a single return 
witn iSero tiitte delay is 


ttr 

O ^ 

Cg(t) * R p ; f(d)f (h~t)dr) for <t < 
t 


» 0 otherv7i&e 


f (7)) becomes zero for ?? <0 and 'll ^ 
*• Tor 0 1 we have 

a 

o ^ 

c„(t) ^ R p / f(n) f (7j-t)d?j 

t 


R.P 


r 


E o 






« R P. 


T 

, ^ e^“°‘ pat(2n_t) 

t_ 
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0 

'T i ^5+at(2Tr.t)l dT) 

0 t 


iT^oTt*’ ^ (2T^~t)3 - sin ] 


E 


X oTt ^ <i't(I ~.t)J Gos(f> t+atT ) 

o ^ o 0 


Similarly for -T <t ;^0 we have 


t+T 


o 


C (t) « R P / f(r,)f {,,»t)d/'l 


t+T 


= R P 




dr? 


O 


t+X 


T 


/ cos t- at(2r?-t)) dr? 


t 


S 

jf M4 ^ 

2T 


T^r [ ‘ ain {C4 t h at(2T +t)J-sin(w t c t^) i 
^ o O O J 


“ ""x^atT^ ^ aos{a)^t -f ^tT^) 

*^0 

A 2 VARIANCE OR THE NOISS PAPT Or THE OUTPUT QP CORRELATOR 

i^cTiVsir 

HU I 

The correlator receiver is a linear time ~ invariant 
iiltor with impulse response h^(t)/ given by equation (5 39) 
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Since the input noise is white v/ath spectral density n/ 2 v;e 
hive 

>» 

E[ I C^j(t) ] = 0^ = I /|h^(t)|^dt (A 1) 


2r 

o 


0 

/ dt 



EN 

"^2 



APPENDIX B 


SIGNAL P/\Ri or IHE OUIPUT OP V® BASED RECEIVER iOR A 

SINGLE iUTORN 


Signal part W (t) of the output of the WD based 
s 

receiver foL a single return with zero time delay is given 
by equation (6 29) 


W„(t) s R P 


r^Tt/2 

/ 

t/2 


j{co +2a(x-t/2)} 


t 


dx 


« 0 otherwise 


for -2T^ < t <2r^ 


(B 1) 


For the sake of conpubation of the above integration the interval 

<t < 2T would be divided into smaller sub-intervals 
o — o 

Gaso 1 

-2T < t 1 -T 

o “ — o 

Since Wv, 6)) is nonzero only for 0 < t < T^, the limits of 

H/ n 

the integration in equation (B 1) become zero to (T^+t/2) 
Substituting the value of ^ equation (5 7) and 

notiag that T^'it/2 < T^/2/ vre get 



R P 


T ht/2 

Jt“o+2Cv(x-t/2)} t 

2T atx ~ e dx 

n 0 


2jgj 

r rt sin(2c.tx:) cos{« tr? (x-t/2)tJ dx 
0 ^ ° 

ht/2 

57 tt ^ { sin(w^t+4o.tx-ct^)-,cin(-o)^tH-<: t^)} dx 

0 

X^Tt ^ ^ (w^t+2atT^_ut^) + X t cos ( 

2 

«cos(-.0' t:'i2c.ti +at )} ] 
o o •' 

j/W [ (Jl' + o) sin((' t+2cuta { sin(-a) t^■v^.tT ) 

rCt.'- 02 O O 2Cfct o o 

0 

sin{a.tT^+ Cut^) } ] 


E(2t +t) 

t ^ _ 

o 


sin(2atT -<i,t^)+ — ^-y-y sin(atT ) 
o 20? a'^t ® 

o 


ein(0'trQ'K‘‘t^) } cos 
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L(2i +t) 2 - 

+{ 2*r~ut! — c:os(?ctT ^ u.t ) «• Gos(cttT )sin(ct.T +o-t )j 

o 20? Ci- t ° 

o 

, ain(6'^t) 


Case 2 

, t < 0 

o S. “* 


Here bho limibs of the integration in equation (3 1) become 
arro to (T^ > t/2) as in the previous case. But since in this 
interval (T^')t/2) >, we get 


0 ? /2 
o 

Wg(t) = R.P / 

0 


4Ex sin(2atx 

2t o.tx 
o 


*"1 ( x — 1 / 2 ) } 

e ° 


t 


dx 


+ R P, 


r^-it/2 


a /2 
o 


4D (1 ~x)sin 2(it{T -x) J (w +2. ( c-.t/2) } 
2T^at(T^«30 


T^/2 


E 


/ {igin(o^^tH4c».tx-<xt )+sin(-/o^t Kvt )) dx 


0 
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T +t/2 

h 9 o 

^ / { sm (o)^t-^2ctx-ut ) sin 2Liut.j.^--.4cLtx <t't )} dx 

1^/2 

+ ^ sin(^)^tKi,t^) ] T Ir t slii(0)^t+2atT^^^t^) 

1 2 2 
H I C0S(-^Qt-2c<.tT^^t )-cos(-t>^tK-L ) }] 

“ gf ’cTt sin(at^) sinij^^t+2atl^} +T^ sin(a.tT^) 

cos (t-' ti-u-tl ~ot^) + sin(w t 2a. tT -l t^) ] 

O 0 2 o o -* 

~{ oos (ctT ~at^)sxn(v-tj?^) 2f~G, sin(2atT^^ t^) 

o 

^ .* -S^-y sin(ot^)sin(2atT^)} cos («^t) 


+{ 


E 


2J? u. 

O 


Gos(2a.ta?^.-u.t^)-. “'■^'-"T“2 sin(at'^)aos(2<atT > 

o 


E 

Wet. 


Sin (k tT^-wt^sinCutT^)) sin (w^t) 
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Case 3 


0 < t < T 

- o 


Hare the or integrution become t/2 to T Also t/2< T /2 

o o 




^^g(t) « R P 


t sl n(2u^ J{o)^+2u( <-t/2)} t 

2T o.tx " ® 

t/2 o 


dx 


+R i> 


^ 4L(T «x) sin 2ct(x -x) J {c.)^+2a(x-t/2)} t 

f vrL s ^ 


T^/2 


2T StO? _x)' 
o o 


dx 


E ? 2 

^ { sxn (o.^ t.+4'u.'Lx^ lat ) sin(^co t <ut )} dx 

t/2 ^ 


o 


T 

O 

L 2 9 

+ V''c’+-'” ^ ^ sin((») t+2atT «Cv.t )tsin(-(i> t+2c-tT -lo-tx+ct )} dx 

O O O O 

X/2 


o 


T 1 

(-^•.) sin(^^tKut^)] + yIt'C 5^ sin(w^t+2u.tT^-Ut^) 


1 2 2 
•' r '- r { COs(-.(j) tKvt -2CutT )-cos(-0' t+at ) }] 
4 utl O O O 
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” 0? L t t t )~COs(-<j) t.+wt^)j+ (-'2^) Sin (-o> t-K t.^) 

o V ^ O 

T 

+ ^ sin ((>3^t+2utT^-ut^) j 


Erl 2 7 

'^’ 1 ) Sin(<-W^t) + (—^■■ 2 ^) sin(— ) 

T 

H Sin ( «^t+2v tT^-. t ) ] 


E (T 


+ {“4 aos(2atT -ct^) 
U*-*C o 


sin. c t‘ 

2r^ c V 


E (i ►-t) 2 

-’ 2 '*|,’^£-"^ cJOs{u.t ) } sin 
■^o 


Case 4 

Here the limit of the integration is t/2 to and t/2 > T^/2 

^ 4B(T -x) -j {OJ t2u.(..^-t/2) } t 

» W_(t) « R P S 2ut(l^-x)} e dx 

a 4*0 

t/2 
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T 

o 

^ {sin(0^t-l-2 tT^-,ut^) sin(-W^t-J-2 t.x^-4 tx-r t^) } dx 

t/2 

r^JlT (WQt-r2ctl^„ th ] 

. [cos(-i^'^t~2v.tT^-K t^)-cos(-tJ^t+2 tr^~3 t^) ] 


[ { r^r-t) s±n.(o^t+2 ti^-v^-t^) l-2sin(-.«^t-. t^) sin(2\. tr^-2 -fc^)] 


E(r ~t) 


r wt 
o 


sin (2 tT^-Lt ) 


L 2 ? 

—•■ Sin (2 tT ~ t-t )sin( fc ) 1 


C06( to 

O 


B (T^-t) 
T Lt 


C0s(2 tT t^)- -y- sin(2i. tT --t^) GOs(>- 1^) } 

o 


sxn{ 

o 



APPEi'IDIX c 


VARIAilGD Ot TUB NOI:,E P/aRj? 0"’ THE OUTPUT OP WD B^ED RECEIVER 

As hds been discussed in section 5 4^ the variance 
of the outout noise in the based racsiver can be computed 
by carrying out the integration of equation (5 34) The range 
R Oi the permitted value's of Tj can be found by dividing the 
intervals x =* 0 to and y = 0 to T^ into aopropxiate smaller 
subin tcrvals Thus the double integiation of equation (5,34) 

can be split into a number of double integrations^ each of 
which will have range R defined in terms of x and y 
Computation of these smaller intervals is straightforward 
since the integrand of equation (5 34) is a just a constant 

The orocedure outlined above is quite lengthy and 
hence will not be included here although it has been verified 
that this leads to the value of variance given in equation 
(5 35) Ins Lead an alternate derivation based as the inpuLse 
responje concept is given below 

As has been mentioned in Soction 5 5, the WD 
bSfBfi® receiver can be replaced by an equivalent linear time™, 
invariant filter v^ith impulse response h (t), given by equation 

rV 
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(5 38) Hence using the formula given by equation (A 1) 


we get 


B [W^(t) f / ^ at 


I (2T^+t j^dt -f ||> 


S T dt 
o 


+~v?r / <at 


■.W [4T^ -- 


p3 OT^ 

^ + ~^3~> 


3 i 3 3,3 


^ + 2T? - 8T^ h ' "^o *’* 3^ ■* ^ 


5T 

Ni- r 1 - 2 . I 

” 2T~'‘ 3 -* 

O 


5HDT 



APPcNDIX d 


IMPULSE HDSPOUoE OP m EQUIVALEIJT TILIER FOR 

the wd based HECDiysR 


The iirpuxse response h (t) of an ecuivalant 

w 

XiZior wnich can replace the WD based receiv-r is the out- 
pu,t of tho V® based roceivur for an input signal r(t} 


Usinq ecfUition (5 2o) we get 
ti-i 

O 1 00 

h^(t) « / r fiv -1 {(iV,^2o( ,„t) } p 


f^fst p»»f«00 


do dr/ 


t+T 


° ^ 2jn {w j-2Cv(d^t) } 

/ 2 £^{ 2 me ^ 


dr) 


t+l' 


« / 2. £ (2n-t) e 


2 jr) h 2 a(o~t) } 


o 


dn 


Now oan be aornputed for different ranges of t 
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C jO 1. 


t < 

AS II goQS from t to (t+T^), goes from t to (1J 2r^) 

Since both t and (tt2T^) are neyitive values we get a zero 
intagrition value, 

c 1 j© 2 


< t 1 

Q ^ 0 

Here the limits of ?7, which vail keep (2T/-t) in the range aero 
to are t/2 ind (ti-T^) 

.M r9n ^j^(2n-t)^ 2jr){(o +2u(?i-.t) } 

j«^(2h-t) ^ e ® dn 



t+T 


bjt) 


o 

S 2 


rr 

T^o 
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Case 3 


-,i < n < 0 

O — 


T 

ticre limit-s of should be t/2 and 


Using equation (D 1) ^/o got 

(t^rt)/2 

h^(t) « / 


^ 


■A 

„ -J&t^ 

T e e 

0 

Case 4 


0 t , 



T 

lleio the required limits on 0 are t and (— 

UoJng equation (D 1) 

we yet 

h^(b) « 

0 ]£ 

/ 2 e 

t >J 0 

=iV 

-jat^ 

. T -t e e 

0 

V ^0 

Case 5 

.4. tt, 


0 

< t 


e 


-jat 


As in oasa X here also there no range of H from t to (t+T^,) 
which will moho f*(2ivt) nonsero Benoe h„(t) =0 



APPENDIX E 


ENVELOPES or JIE OUTPUIS OP THE CORRELATO-l 
RvCi^lVER AND TIE V® BASED RE CL TVER FOR Ol-d 
RETURN : SIMULATION RESULTS 
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APPDNDIK F 


^w^LOPus or ri-E oaieurs or ihc corr!i.l.ator 
RECEIVER AND TIC WD BASED REC«1V-,R POl TWO 


RETURNS I SIMULATION RESULTS 
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